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IXKXAF  FROCRAMira 

by  Va.  Orchard -Bays  (Rand  Corp.) 
Hal  Judd  (I.B.M. ) 

Leo  la  Cutler  (Rand  Corp.) 

nraopqcTioK 

ft*  linear  programing  lyitn  vhlch  fell eve  la  the  modified  simplex  pro- 

'\ju  jUusC* 

cedure  vlth  the  product  fora  of  Inverse^  It  la  designed  to  solve  the  class  1- 
cal  linear  lnequa''  It  lea  problem  and  aort  of  Its  variat  ions  on  the  IBM  704.  / 

jftgnai  w&nxggm 

The  basic  Machine  needed  for  this  systm  Is  the  IBM  704  vlth  4  logical 
dnnm,  4  tapes  and  4096  wards  of  Magnetic  core  storage.  One  additional  tape 
■ay  be  used  optionally  for  storing  output  results.  One  tape  aay  be  used  for 
input  to  the  data  as  seat)  ly  vhlch  Is  run  prior  to  using  the  system.  Hence  the 
input  is  from  oards  or  tape  and  the  output  is  the  on-line  printer  or  off-line 
printer.  The  card  punch  Is  used  by  the  data  assembly  and  also  for  punching 
restart  Information. 


Mormon 

The  quantities  displayed  In  a  Matrix  require  tvo  indices,  one.  for  rev 
and  one  for  colon*.  Here  a  row  vector  vlll  be  denoted  by  a  general  subscript, 
one  of  lte  elenents  by  a  specific  subscript,  as* 

Cj  Is  a  row  vector 

c^  Is  the  el  went  of  Index  5. 

A  colon  vector  vlll  be  denoted  by  a  general  superscript,  one  of  Its 


elements  by  a  specific  superscript,  as: 


bJ  is  a  colon  vector 


b*'  Is  the  element  of  index  5. 


A  natrix  vlll  be  denoted  by  both  a  subscript  and  a  superscript,  as: 
a*  Is  a  natrix 

a^  It  the  row  consisting  of  the  element  of  ludex  5  in  sach  colon. 
J 

I 
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a?  is  the  colunn  consisting  of  the  element  of  index  4  la  each  rev. 
fc':;  le  the  element  la  roe  6  and  colimn  2. 


The  ranges  of  values  will  be  specified  when  the  index  is  first  Mentioned. 

Hots  that  a1  Is  the  transpose  of  a^.  The  transpose  of  a  na  ■  lx  Bust  be 
denoted  by  defining  a  new  letter,  e.g.,  b*  =»  af  defines  bj  as  the  transpose 

For  a  given  eat r lx  of  coefficients  a*  and  a  colusn  of  constants  b1,  then, 
a  system  of  n  simultaneous  linear  equations  in  n  unknowns  x?  will  be  written 


a*  x^  *  b‘ 


(l  =  1,  •••,  ■!  J  *  1,  «..,  n) 


instead 


of  f 

J-T 


xJ .  b1 


for  1*1,  ...,  n.  This 


tlon  convention  is  used 


whenever  a  column  is  Multiplied  on  the  left  by  a  row.  A  Matrix  nay  be  Multiplied 
on  the  left  by  a  row  vector  to  produce  a  row  vector: 

Ci  "  dj  * 

A  matrix  nay  be  Multiplied  on  the  left  by  a  row  vector  and  on  the  right  by  a 

column  vector  to  produce  a  scalar: 

1  J 

ci  ‘j  *  -  *  • 

Finally,  of  course,  a  Matrix  nay  be  Multiplied  on  either  side  to  produce  another 

■atrlx  by  Multiplying  by  a  satrlx  with  the  proper  dimensions: 

1  VJ  1 
•j  \  -  °k 

.h  .1  h 
4t  ’  eJ  1 

Uppercase  letters  (except  T)  will  be  used  for  sets.  The  use  of  letters 
for  Indices  is  as  follows: 


h,  1,  J,  k  for  general  Indices 

I,  n,  n  for  limits 

p,  <b  T>  ®  for  specific  indices. 

When  it  is  desired  to  Index  a  quantity  over  time,  l.e.  iterations,  the 
index  will  be  encloeed  in  parentheses,  e.g.,  is  produced  from  after 

t-1  transformations.  The  capital  letter  T  will  be  used  as  the  current  limit  for 
t,  that  is,  T  Is  the  current  iteration  mmiber  and  t  =  1,  2,  ...,  T. 


mam  <*  nctim 
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V©  define  the  standard  linear  programing  p rob  1m  In  tvo  for**.  Thle  Is 
because  It  le  Mandatary  to  start  vlth  the  identity  Matrix  as  the  initial  basis 
in  all  clrcuMtaneee.  Certain  deviations  frcn  these  standard  fores  are  dis¬ 
cussed  la  the  detailed  vrite-up. 

BtAMPABP  FCRM  1 

Given:  A  rev  of  coefficients  01  a  linear  foam  to  be  opt  ini  ted 

•5  <> 1 . *> 

A  BBtrlx  of  restraint  coefficients 

aj  (i  *  1,  s;  J  -  1,  ...»  I) 

A  col  urn  of  constants 

b1  (i  -  1,  n) 

To  find!  A  col  vast  of  values  far  x^  ^  0  (j  «  1,  . t)  such  that  the 
variable  x°  Is  Marini  ted  subject  to 

(1.1)  *  +  aj  ^  “  o  (j  3  I#  0 

(1.2)  a*  xJ  £b*  (1  -  1,  n) 

Vote  ti*t  (1.1)  1s  perfectly  general  since  the  sun  a®  x**  nay  be  Minimized  or 
■ax In  1  zed  Merely  by  changing  the  signs  of  the  a®.  To  convert  (1.2)  to  equali¬ 
ties,  ve  define  the  variables  x*+*  si  0  (l»  1,  ...,■).  Then  (1.2)  beecnee 
(1.2*)  *J  *  x'*1  -  b1  . 

If  8^  (i  -  0,  1,  ...»  a;  h  =  0,  1,  ■)  is  the  identity  natrix  of  order 

n+l  (essentially  the  Krooecker  delta),  then  ve  can  define 

11  o 

a  -  b  and  b  -  0 
o  o 

Vh  "  *h  for  h  “  lf  •*'  ■ 

(thu*  defining  J  -  0  on d  )  -  1+1,  Im  -  b)  and  replace  both  (1.1)  and 
(1.2')  by 

(1.5)  (l  -  o,  1,  ....  ■;  J  =  0,  l,  n). 

If  the  restraint  equations  can  be  put  in  the  for*  (1.2' )»  then  sene  of  the  b* 
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are  permitted  to  be  negative,  if  neceaaary.  Bovever  the  number  of  nsfatlve  b 
should  be  Mall  to  prevent  as  excessive  maker  of  iterations. 

If  the  restraints  cannot  be  specified  exactly  in  the  fom  (1.2)  without 
an  excessive  maker  of  negative  b1,  then  the  problsa  should  be  put  in  standard 
fora  2.  It  should  be  noted  that  the  slack  vectors  Introduced  above  are 

legit  lasts.  The  Halt  n  will  be  used  consistently  for  the  maker  of  legit  lasts 
vectors  in  the  system  (besides  a*  which  is  an  operational  device.) 

BTAflDAHD  PORK  2 

Olvent  A  row  of  coefficients  of  a  linear  fora  to  be  optimised 
*j  (J  e  1#  ...»  ») 

A  aatrix  of  restraint  coeff  ici^^ts 

aj  (1  -  2,  3#  a;  J  -  1,  2,  ...»  n) 

A  coliam  of  constants 

t'jo  (1.2,  5,  a) 

To  find :  A  coltau  of  toIum  for  r  io  ( J  •  1,  . . . ,  n)  auch  that  th. 
variable  x°  Is  aaxlmlzed  subject  to 

(2.1)  *°  +  »°  t3  -  0  (J  -  1,  ...,  n) 

(2.2)  .J  xJ  -  b1  (1  -  2,  ...,  ■)  . 

Ve  now  add  to  the  system,  artificially,  the  Identity  aatrix  (except  6* 
which  Is  always  there)  and  auxiliary  variables  x  (k  «  1,  ...,  a)  and  con¬ 
struct  an  auxiliary  fora  in  which  the  variable  x0*1  Is  to  be  maximized  first, 
l.e.  before  maximizing  x°.  This  process  is  called  Phase  I. 

The  purpose  of  Phase  I  is  to  eliminate  the  artificially  added  coltame 

from  the  system  or  at  least  to  aake  sure  that  the  corresponding  variables 
# 

become  tero.  Bence  ve  put  a  vdlght  of  unity  on  each  one  in  the  rw  of  Index 
1,  vhich  will  become  the  auxiliary  form.  Let  =■  and  for  k  ■  2,  ...,K) 

*  'At  least  one  artificial  coliam,  usually  the  coliam,  must  remain  in  the 
solution  at  zero  level  even  while  maximizing  the  variable  x°. 
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l*t  a^>jc  -*  +  Th©  variables  x11^  are  to  be  non-negative  except  for  a?1** 

which  it  defined  by 

(2.3)  x**1  + 

k»cF 

Clearly  xnf^^  0  and  if  x11**  =0,  than  all  x04^  „  Vhen  (and  if)  thie  con¬ 
dition  it  attained,  the  artificial  variables  a re  maintained  at.  zero,  including 
# 

x  ,  by  considering  (2.3)  at  a  restraint  equation  while  Maximizing  x°  (Photo  II.) 
If  xml  cannot  be  driven  to  zero,  there  it  no  solution  to  (2»2)9x‘^tO.  This  con¬ 
dition  it  called  Terminal  1. 

The  problem  at  this  point  can  be  displayed  in  the  following  augmented  form 
where  the  x^  are  thown  above  the  matrix  of  detached  coefficients. 


Since  the  identity  matrix  rtill  does  not  appear  in  (2.4),  ve  must  make  a  simple 
preliminary  transformation.  First  note,  however,  that  setting 
x  *  0  for  J  ■  0,  1,  ...,  n 

n+1  yS  1  n+i  .  i  .  ,  «  _  _ 

x  »_2»b  and  x  =»  b  for  i  »  2,  3,  ...»  m 

i=£ 

n-f  1 . 

provides  an  initial  solution  in  which  all  variables  are  non-negative  except  x 
This  solution  will  remain  valid  If  we  subtract  all  equations  of  index  2  through 


m  from  the  auxiliary  optimizing  form,  row  1.  Then  defining  anew 


1 

J 
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(2.6) 


for  J  *  1,  2,  ...,  n; 


0  for  k 


2#  5>  •••»  * 


the  rov  1  equation  takes  the  fora 

(2.7)  xJ  *  I0*1  »  b1  (J  -  0,  1,  ....  n). 

furthermore,  the  colusns  s*  ,  a^,  a^2  ,  a*^  nor  fora  Bj^  . 

The  data  assembly  program  computes  aj  and  b1  automatically.  If  some,  but 
not  all,  of  the  coluns  of  B^  occur  In  aj  (j£n),  then  these  colimns  should  be 
indicated  as  being  in  the  initial  basis.  The  data  assembly  vill  then  omit  the 
corresponding  rove  in  the  suae  (2.6).  This  is  equiralert  to  adding  artificial 
columns  a*+k  only  for  those  colusms  of  B^  vhlch  are  missing. 

The  rectors  a*  ,  a*^  are  never  entered  with  the  data.  They  are  implicit 
in  the  code  and-,  once  eliminated,  cease  to  exist.  The  basis  headings,  vhlch 
are  the  J-lndlces  of  the  columns  a^  in  position  h  *  0,  1,  ...,  a  of  the  basis, 
are  left  zero  for  all  artificial  rectors  since  the  position  h  identifies  them 
sufficiently,  l.e.  they  are  merer  moved  around  in  the  basis.  Legitimate  col¬ 
umns  of  B^  nay  go  out  of  the  basis  and  come  back  in  out  at  position.  Bence 
they  must  hare  names,  as  In  Standard  Fora  1. 


Note:  The  rov  indices  1  of  a*  must  hare  the  numerical  raluen  0,  1,  ...,  m, 
but  the  column  indices  J  =  1,  ...,  n  are  used  above  only  for  expository 
purposes.  These  columns  can  be  identified  by  any  n  distinct  symbols,  of  fire 
or  less  Hollerith  characters  each,  which  suit  the  fancy  of  the  f emulator  of 
the  problem.  One  convention  has  been  adopted:  legitimate  columns  of  B^ 
are  denoted  by  DP001,  UP002,  ...,  for  "Obit  Positlre”  vector.  This 
convention  is  essential  only  to  the  re-inversion  code. 
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Yh«  data  asseekly  far  the  linear  programing  system  la  doalgnad  to  be 
executed  Independently  of  the  win  code.  The  opt  lone  for  loading  data  from 
tape  or  from  the  card  reader,  to  load  uatrlx  elenente  vhlch  are  punched  In 
etandard  font  or  In  a  fixed  field,  or  for  punching  out  the  natrlx  on  binary 
card*  are  not  controlled  by  eenae  evltchee.  Instead,  each  bit  of  a  vord  In 
storage,  the  81NBXB  vord,  la  ueed  to  control  these  functions.  In  particular, 
the  flrat  three  bite  correspond  to  the  three  options  noted  above.  A  zero 
vould  correspond  to  the  switch  being  up,  and  a  one  to  the  switch  being  down. 

A  binary  card  suet  be  Inserted  1  need  lately  following  the  transfer  card  In  the 
data  asembly  deck  with  the  8XNBB  vord  punched  In  the  9 'a  left  row.  If  the 
data  le  loaded  from  the  card  reader,  a  one  It  punched  In  position  S  of  the 
SmS  word,  l.e.  a  9-punch  In  colvnm  1  of  the  card. 


The  xmm  of  the  problem  or  any  suitable  Identification  le  punched  In 
eoliams  2  through  72  in  Hollerith  characters  for  the  heading  card.  This  Infor¬ 
mation  Is  punched  In  a  binary  card  along  with  other  data  for  use  by  the  main 
code  vhlch  prints  the  informat  Ion,  as  origlna^punched;  as  the  flret  line  on  all 

of  the  output  list  Inge. 

# 

raOBUPt  PARAKCTR8 

The  essential  parameters  needed  sure  punched  in  the  parameter  card  In  nor¬ 
malized  font  far  ease  In  key  punching,  e.g.,  a  equal  to  77  vould  be  punched 
In  eoliams  1  and  2  while  n  equal  to  177  vould  be  punched  In  columns  1,  2  and  5. 
If  the  makers  are  not  In  normalized  farm  on  the  card,  leading  zeros  must  be 
punched.  The  parameter  card  ie  as  follows: 

cole  1-10  a  (maker  of  restraint  equations)  rnZX  7n  »  255 

* '  11-20  z  (number  of  Phase  I*s  desired) 

* '  21-50  q  (total  maker  of  Phases) 

"  31  -  40  £  (index  of  sia  row) 

"  41-50  t  (maker  of  preliminary  transformations) 
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It  the  sub  row  is  computed  prior  to  assembly  and  loaded  with  the  ether 
elements,  T  should  be  left  blank.  This  Is  necessary  if  sore  than  one  Phase  I 
le  desired  since  there  will  be  as  aany  auxiliary  ferns  as  there  are  Phase  Z's. 

If  there  is  no  sw  rw,  leave  £  blank,  not  zero. 

BASIS  HIAPIWOS 

The  basis  headings  are  punched  7  to  a  cord  using  as  aany  cards  as  needed. 

Ten  colxssns  are  allotted  for  each  heading  ( j)  with  Its  colmn  Index  (l).  The 

colvasn  index  Bust  be  in  nomall  zed  form  or  leading  zeros  Bust  be  punched.  The 

cards  are  punched  as  follows: 

cols.l  -  3  1  (coluan  index  far  basis  heading) 

col.  4  blank  or  Blnus  (a  sinus  sign  will  cause  the  sain  cade 

to  treat  colvsm  1  of  the  initial  basis  as  an  arti¬ 
ficial  unit  vector,  and  then  to  bring  the  specified 
rector  j  into  the  basis  arbitrarily.  The  number  of 
negative  basis  headings  Bust  agree  vlth  paraaeter  t.) 
cols. 5  -  9  J  (5  Hollerith  characters  to  denote  ccluas  1  of  the 

identity  aatrix.  These  unit  vectors  Bust  be  denoted 
t.  UPOOl,  UP002,  etc.  for  the  re  inversion  code.) 


col. 10 

blank 

cols.  11-20 

oa  cols.  1-10 

cols.  6l  -  70 

as  cols.  1-10 

If  one  of  the  basis  headings  Is  punched  erroneously,  it  can  be  overwritten 
without  repunching  the  entire  card  by  succeeding  it  vlth  a  correctly  punched 
version  of  the  heeding  in  error.  A  heading  of  zero  will  be  assigned  far  all 
i's  not  specified. 

RIGHT  HAND  SIDE 

The  right  hand  side  vector  b*  is  preceded  by  a  card  with  FIRST  B  punched 
in  Hollerith  characters  In  coluans  1-7,  followed  by  elenent  cards  for  non¬ 
zero  entries  as  follows: 

cols.  7-9  1  (row  index  in  nornalized  form  or  with  leading  zeros) 


p 


.? 


col*  10  sign  of  elcanst  (blank  is  interp* eted  as  positive) 
col •  11“*  The  integral  part  of  the  nuariber  is  punched  in  col.  11  on 
depending  on  the  mmfcer  of  columns  needed  (8191  is  the 
maximum  Integer  allowed).  The  integer  is  separated  froa 
the  fraction  by  a  blank  coltan.  Then,  the  fractional  pert 
is  punched  using  as  nany  coliams  as  needed  (a  maximum  of 
seven  coltmms  nay  be  used).  The  regaining  columns  on  the 
card  are  left  blank.  For  example,  .0 k25  vould  have  0b25 
punched  in  cols.  12*15.  Also,  10.0  vould  have  10  punched 
in  cols.  11*12  only. 

An  alternative  method  has  been  provided  to  load  data  which  Is  punched  in 

a 

fixed  eoltmns  as  described  below.  A  one  is  punched  in  position  1  of  the  SENSES 

word  if  data  is  to  be  loaded  in  this  fora. 

coir.  11-lU  Integral  part  of  mnber 

col.  15  blank 

cols.  16-21  fractional  part  of  nuaber 

If  a  supplementary  right  hand  side  c*  le  to  be  loaded,  a  card  vlth  NOT  B 
punched  in  col  tame  1*6  precedes  the  element  cards  which  are  in  the  same  fora 
as  ths  elanents  of  to1. 

m 

A  FIRST  B  card  must  always  be  used  even  if,  far  some  reason,  no  b1  are 
entered,  e.g.  when  assembling  a  supplementary  right  side  or  a  new  matrix. 

MggtADjT  MATRIX 

A  card  with  MATRIX  punched  in  colusns  1*6  precedes  the  element  cards  of 

the  natrix  a1.  If  no  natrlx  is  to  be  loaded,  this  step  is  bypassed  by  loading 
J 

a  card  with  90  MATRIX  punched  in  colimms  1-9.  The  non-zero  elements  of  the 

matrix  a*  are  loaded  one  element  per  card  as  follows} 

col.  1  always  blank  (corresponds  to  sign  column  in  basis 

heading  cards,  but  must  be  plus  here.) 
cols.  2-6  J  (a  five  Hollerith  character  symbol  to  identify 

the  column).  All  of  the  elements  of  each  column  vec¬ 
tor  must  be  in  succession. 


r-oe* 


P.10 


cole.  7— >  row  index  and  element  punched  In  the  ease  fora  at 

i  • 

eloents  of  b  . 

A  card  vlth  SCP  punched  In  coliame  1-3  will  terminate  loading  of  the 
matrix.  A  card  vlth  ECR  In  colunel-3  will  cauee  the  data  assmribly  to  write 
the  preceding  Teeters  on  the  tape  and  start  a  new  record  vlth  the  succeeding 
vectors.  If  position  2  of  the  SXT6B  word  Is  a  one  the  data  assembly  punches 
the  matrix  In  binary.  When  It  Is  loaded  by  the  main  code  It  is  possible  to 
load  part,  but  not  all,  of  the  vectors  If  they  have  boon  separated  by  the  SOR 

breakpoint  (see  Deviations  from  the  standard  farms.)  Also,  a  card  vlth 

# 

CURTAIN  punched  In  coltmms  1-7  may  be  loaded  between  any  tvo  vectors  (see 
Deviations  from  the  standard  farms.) 

If  a  sum  row  Is  used  for  Phase  I,  then  the  sum  of  each  colimm  must  be  less 
than  8192.  This  does  not  apply  to  the  b*. 


* 
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DPAIIP  VRITI-U? 

The  codes  vhoee  use  ere  described  herein  ere  designed  to  provide  e  com~ 
plete  system  for  solving  linear  programming  and  related  p  rob  leas  where  no 
special  assumption*  are  node  regarding  the  structure  of  the  model.  Since 
there  are  so  many  ratifications  to  such  a  general  c caput er  program,  even 

for  one  type  of  proble*,  it  Is  desmed  necessary  to  Include  a  considerable 
aaouxrt  of  explanatory  naterlal.  In  vhat  follows,  a  general  familiarity  with 
the  sljgplex  method  will  be  assumed.  For  example,  no  proofs  will  be  given 
for  the  standard  theorems  concerning  basic  solutions  or  the  simplex  criterion. 
Bower er,  toe  algorithm  using  the  product  for*  of  inverse  and  certain  special 
features  built  Into  the  codes  will  be  developed  In  detail. 

ftie  subject  matter  will  be  divided  Into  the  following  sections. 

I  Advantages  of  the  product  far*  af  inverse. 

II  Reducing  the  nxmiuer  of  terms  in  the  product.  (Pe-lnverslon) 

III  Deviations  f'-om  the  standard  farms. 

IV  Farms  of  arithmetic . 

V  The  basis 

VI  The  Inverse  of  the  basis.  (Product  Form) 

VII  The  pricing  operation. 

VIII  Choosing  the  basis  coluan  to  be  replaced.  (Degeneracy) 

IX  Summary  of  cyclic  operations.  (One  iteration) 

X  The  composite  algorithm.  (Infeasible  solutions) 

XI  Parametric  programming.  (Varying  the  right  hand  side) 

XII  Multiple  phaeee.  (Varying  the  optimising  form) 

I  -  ADVANTAGE  OP  THE  ITCDtET  FORM  OF  IKVPfil 

In  order  to  appreciate  the  advantages  of  this  method,  it  Is  necessary  to 
keep  in  mind  the  following  facte  concerning  the  matrix  of  coefficients  of  a 
typical  linear  programing  p rob lms  and  the  requirements  of  the  simplex  method. 
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(A)  The  matrix  Is  usually  Tory  sparse,  the  percentage  of  non- zero  ele¬ 
ments  being  typically  5  to  1?  per  cent  of  the  total  number. 

(B)  The  given  data  Is  almost  Invariably  expressed  vlth  a  very  fev  signi¬ 
ficant  figures  and  Is  easily  scaled  to  fit  in  a  restricted  range  of 
fixed-point  numbers. 

(C)  Scam  prob learn  have  a  large  mxZber  of  variables  although  the  xmafcer 
of  restraints  Is  not  excessive.  The  ratio  of  nxzfcer  of  variables  to 
number  o'*  restraints  Is  often  betveen  5  to  1  and  10  to  1  and  actual 
problems  have  been  run  vlth  ratios  as  high  as  50  to  1. 

(D)  If  m  Is  the  number  of  restraint  equations,  then  the  essential  manipu¬ 
lations  of  one  cycle  of  the  simplex  method  involve  a  change  of  basis 
in  Euclidean  m-space,  where  two  successive  bases  differ  only  by  one 
column  vector.  The  rules  of  elimination  in  effecting  this  transfor¬ 
mation  are  the  same  no  matter  vhat  particular  algorithm  is  used.  Al¬ 
though  such  a  transformation  is  easily  represented,  its  total  effect 
on  the  mamrlcal  representation  of  a  set  of  vectors,  l.e.  a  matrix, 
may  be  extensive.  Furthermore,  the  validity  of  each  cycle  performed 
depends  on  the  accuracy  of  all  preceding  cycles. 

Vlth  these  observations  In  mind,  tho  advantages  of  the  present  method  can 
be  stated  as  follows : 

(1)  Since  the  original  data  natrlx  is  not  transformed  from  iteration  to 
Iteration,  It  Is  clear  from  (A)  and  (B)  that  an  elaborate  organiza¬ 
tion  of  the  data  can  be  performed  at  the  outset,  so  that  it  Is  In 
the  most  convenient  and  compact  form  for  use  throughout  the  probltm. 
Since  the  original  matrix  is  referred  to  only  once  during  the  cycle, 
it  can  be  stored  on  tape,  out  of  the  way,  and  its  compact  form  mini¬ 
mizes  transfer  time  as  well  as  computing  time.' 

(2)  It  is  clear  frcm  (C)  that  it  would  be  very  expensive  in  some  problems 
to  transform  the  vholo  natrlx  an  each  cycle  since,  even  if  the  ratio 


of  variables  to  constraints  Is,  say,  only  as  much  as  5  to  1,  it  is 
•till  rory  unusual  for  the  process  to  take  rjm  iterations.  In  other 
vords.  Many  rectors  (activities)  are  never  selected  at  all  hen; a 

there  is  no  use  doing  anything  to  them. 

(3)  If  a  full  inverse  re  re  Maintained,  it  vould  be  necessary  to  store  it 

by  revs  nr  else  linear  combinations  of  revs  would  not  be  readily  avail¬ 
able.  fven  then,  one  would  requiro  storage  for  three  vectors  slmul- 
taneously  in  transforming  a  column.  In  changing  basis,  the  transfor¬ 
mation  vector  is  a  col  tom  and  hence  the  Modification  and  re-recording 
2 

of  the  m  elements  would  be  awkward  and  tljee-ccnsumlng,  especially  If 
zeros  are  deleted.  On  the  other  hand,  the  product  form  requires  the 
recording  of  only  one  additional  column  (plus  an  index)  on  each  iter¬ 
ation.  These  columns  will  alaoet  surely  contain  many  r.eros  and  nay 
be  condensed  since  they  are  always  applied  coluMn-vise  and  recursively 
to  a  single  vector. 

(U)  The  product  form  is  extremely  amenable  to  modifications  of  the  method 
or  for  generating  additional  side  information,  since  the  inverse  or 
lte  transpose  are  equally  easy  to  apply. 

ii  -  RiDTEiaa  Tg»  eucbr  or  tbwb  in  the  racDocT 

One  apparent  d  1  sad va stage  of  the  product  form  is  that,  while  only  one  new 
column  is  recorded  each  cycle,  it  is  necessary  to  read  T  coltsmis  to  apply  this 
form  where  T  Is  the  number  of  iterations  already  performed.  Thus,  after  m 
Iterations,  where  a  is  the  number  of  restraints,  one  must  read  mare  information 
than  In  reading  a  full  inverse  (not  taking  condensation  into  ac court . )  This 
is  still  profitable  for  something  error  2m  iterations  as  con  be  seen  as  follors. 
With  an  explicit  inverse,  we  must  read  it  twice  and  write  it  once 

each  iteration,  giving  columns  handled  on  T  iterations. 

With  the  product  form,  we  must  read  t-1  columns  twice  on  iteration  t 
write  one  new  column,  giving 
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£  2(t-l)  ♦  1  *  T(T-l)  ♦  T  «  colxmns  handled  on  T  Iterations. 
t«l 

totting  T~  >  WT  gives  T  -  >.  However  something  aust  be  allowed 
far  the  fact  that  lets  arithmetic  accompanies  the  writing  than  the  reading. 
This  disadvantage  Is  aare  apparent  than  real,  hoverer.  After,  say,  2a  Itera¬ 
tions,  round-off  error  will  begin  to  be  noticeable  an  large  probleas  no  natter 
which  fora  Is  used.  Many  problems  are  solved  In  2m  iterations  or  less  but,  if 
It  takes  more,  one  can  re- Invert  the  basis  matrix,  at  any  tlae,  producing  not 
aore  than  a  columns  of  Inf oraat Ians.  The  tine  for  this  inversion  Is  auch  less 
and  the  accuracy  of  the  resulting  transformations  is  as  good  or  better  than 
after  the  saas  nualber  of  full  simplex  cycles.  (The  order  of  elimination  for  in¬ 
verting  Is  designed  to  aalxxtaln  accuracy.)  A  special  code  Is  provided  for  this 
purpose.  It  Is  useful  for  solving  any  system  of  linear  equations,  especially 
vhere  several  right  hand  sides  are  used  with  a  sparse  matrix. 

It  Is  helpful,  with  t  form  Of  inversion,  to  consider  the  matrix  of 
coefficients  of  the  restraint  equations  as  a  collection,  or  eet,  of  colun  vec¬ 
tors,  Ignoring  the  fortuitous  orde  tng  given  to  the  variables  In  the  formula¬ 
tion  of  the  model.  Whatever  scrambling  of  columns  that  may  occur  In  the  pro¬ 
cess,  Is  recorded  in  a  list  of  basis  headings  which  accompany  and  Identify  the 
basic  variables  of  a  particular  solution.  This  point  of  view  is  adopted  through¬ 
out  the  present  dlecussion. 

hi  -  DimTioie  n<M  the  standard  tong 

Some  discussion  of  alternate  ways  of  starting  a  problem  is  indicated.  In 
Standard  Tarn  1,  the  Inequalities  were  converted  to  equalities  in  the  usual 
way  by  adding  slack  vectors.  The  matrix  then  contains  the  cosq?lete  Identity 
matrix  (not  necessarily  in  proper  columi  order  in  a*)  and  there  la  no  diffi¬ 
culty.  The  variable  corresponding  to  colu»  1  of  the  identity  matrix  is  set 
equal  to  b*  and  the  name  of  the  variable  is  recorded  in  the  i-th  word  of  the 
list  of  basis  headings.  This  Is  dene  automatically  by  tho  data  assembly  pro- 


gran.  The  problem  is  then  ready  to  start  vlth  this  initial  solution.  Ve  con¬ 
sider  a  Solution"  to  consist  of  the  list  of  basis  headings  and  the  values  of 
the  corresponding  basic  variables. 


In  the  usual  set-up  of  the  simplex  Method ,  It  Is  asouned  that  all  ele- 
■ents  at  the  right  hand  side  are  non- negative  so  os  to  Insure  that  a  starting 
solution  vill  be  feasible,  l.e.  non-negative.  This  restriction  tee  been  re- 
noved  In  the  704  progress  by  Incorporating  a  coenobite  algor  it  hn  which  vill 
reneve  Infeasibilities  as  veil  as  opt  inlee .  (See  section  X  below. )  Thus  it 
le  passible  to  Multiply  through,  by  minus  one,  equations  vhlch  have  negative 

slacks  In  order  to  sake  the  alack  vectors  positive,  even  though  the  starting 
solution  Is  thereby  rendered  infeasible.  However,  it  Is  recommended  that  the 
number  of  Initial  lnfeaeibllitlos  be  kept  saall  to  avoid  an  excoeeive  number 
of  Iterations.  There  le  no  hard  and  fast  rule  about  this  since  different  prob- 
lsas  vill  react  differently  to  the  sane  conditions. 

'If  a  given  Matrix  does  not  contain  all  columns  of  the  identity  Matrix 
legitimately,  l.e.  as  positive  slacks  vlth  taro  coefficients  In  the  optimizing 

l 

fern,  then  ene  can  use  Standard  form  2  and  have  the  code  perform  a  Phase  I. 

In  this  case.  It  le  necessary  that  all  elements  b*  are  non-negative. 

Ixper  lance  has  led  to  the  Incorporation  ef  still  another  device  far  ob¬ 
taining  initial  solutions.  Sometimes  tho  f emulator  of  the  problem  knows  a 
feasible  basis  other  than  the  Identity.  Provision  le  nade  far  Introducing, 
arbitrarily,  any  nxsSber  of  colvam  vectors  Into  the  basis  at  the  outset,  vlth 
the  machine  — the  decision  as  to  vhlch  column  of  the  baeie  each  should 
occupy.  If  the  fomulator  has  misjudged  and  the  specified  columns  produce  a 
singular  matrix,  the  Machine  prints  outs  (and  saves)  the  pertinent  information 
and  stops.  If  the  resulting  Matrix  is  non-llngular  but  the  solution  II  infeas¬ 
ible  (partly  negative),  then  the  casposlte  algorithm  automatically  cute  in  and  works 
toward  feasibility  In  suooeedlx«  iterations.  This  device  ef  arbitrary  trans¬ 
formations  Must  net  be  used  In  Phase  I. 


One  other  provision  of  a  less  obc.lutr  nature  has  bees  node .  Occasionally, 


a  model  la  a  re-work  of  an  older  problem  00  that  acoethli*  la  known  about  lta  be- 
har  lor.  At  other  time  a  the  formula  tor  h&e  certain  Insights  vhlch  he  vould  like 


to  exploit  without  cewrttting  himself  to  absolute  statements  regarding  feasibil¬ 
ity  or  singularity.  It  may  be  possible  to  assemble  the  matrix  columns  In  order 
of  decreasing  likelihood  of  use,  that  is,  vlth  the  moot  likely  candidates  for  en¬ 
try  Into  a  feasible  or  optimal  bat  is  read  in  f  L.  st ,  These  can  be  separated  from 
the  others  by  e  ’burtaln"  vhlch  is  equivalent  to  the  following  instruction:  M  If 
any  candidate  for  entry  into  the  basis  is  available  ahead  of  the  curtain,  use  it) 
otherwise  proceed  to  the  others.  M  Several  such  curtain  narks  nay  be  used.  They 
are  activated  by  a  switch  so  that  their  use  le  under  control  of  the  operator. 

They  have  often  reduced  the  mnfcer  of  Iter  at  ions  required  but  when  used  injudi¬ 


ciously,  l.e.  an  a 


hunch,  nay  have  the  opposite  affect. 


There  are,  of  courts,  other  devices  which  require  no  special  provisions  In 
the  code,  except  perhaps  for  loading.  Far  example,  certain  act ivit lea  (onlnmttf) 


may  to  wltHsald 


the  nubia#  until  optimal  tty  in  dbtaimel  with  the 


of  th#  revised  method  *1  additional  dots  bn  stloj 


«t  any  tine  almply  by  reloading  the  data 


The  product  form  of  lmroreo  In  gaaorated,  mom 
In  Inn  (If),  floating  point  arithmetic.  Certain 


In  »  flmd  1  print  for 


and  cps^Ud  in  da^Mfe* 


then#  regatta  ftto 


tlam. 


*.■*»::* 


put  la  in  flmsd-pnimt  fm  for  mg  mediae,  vlth  6  digits  of  whole  man ff.’lfA 
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We  will  not  bo  so  mich  concerned  with  the  bools  on  itorotion  T  to 


i(T) 


with  lto  inverse  whloh  will  bo  denoted  by  9  ,  that  lo 

(6.1)  ^  £T)  *  ^(T)  %{7)  *  *£  • 

Tho  matrix  novor  exists  explicitly  but  lo  corrlod  u  o  product  of  ele¬ 

mentary  Matrices  which  are  stored  In  a  condensed  for*  consisting  of  at  Most 
one  column  plus  an  Index  and  Identification,  These  columns  are  themselves 


condensed  with  only  non-aero  eleaents  recorded.  They  arc  called  transformation 

i(») 


vectors  (sometimes  elimination  "  equations"  ,  see  (6,5)  )  and  denoted  by  \ 
Since  the  Index  r  Is  Itself  a  function  of  t,  It  should  be  written  r(t).  However 
this  will  not  usually  be  done  since  it  sometimes  appears  as  a  superscript  In  an 
array  already  Indexed  (t-1)  and  confusion  would  result,  ^ 

As  already  indicated,  Oreek  letters  will  be  used  for  a  basis  and  Its 
inverse,  and  also  for  all  (m^l) -order  vectors  expressed  in  terns  of  the  basis. 

s 

Likewise  certain  ratios  and  functions  Involved  in  decisions  concerning  the 
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basis  will  be  denoted  by  Greek  lettere. 


An  explanation  of  the  generation  and  uae  of  the  rector*  le  in  order* 

Since  w*  started  with  for  the  first  iteration,  we  will  hare 

on  hand  to  start  iteration  T.  Suppose  that  coluan  a*^  hae  been  chosen  frcei 
to  replace  coliaui  h  ■  r(T)  in  producing  a  new  baeie  Through- 

e 

out  the  following  discussion,  the  specific  index  r  is  to  be  understood  a*  r(T). 
Superscripts  in  parentheses  will  always  refer  to  the  entire  array  to  which  the 
element  belongs  and  are  not  to  be  understood  as  Modifying  the  open  superscript. 
Let  satisfy  the  equation 


(6.2) 

1 


i(T-l)  h  .  i 
»*h  a.(T)  S(T)  * 


Clearly  CAn  b*  obtained  by 


How 


,,  ...  l(T-l)  h  i 

(6.3)  V  »,(T)  ae(T)  . 

,J(T)  £  .  *(W>  for  b  *  r(T) 


(6.4)  nk'  ’  °h  -  % 
since  only  coluan  h  •  r(T)  changes.  Howerer,  fro*  (6.2)  we  oan  solr*  for  the 


exceptional  coluan. 


(6.5) 

where 

(6.6) 

*/  ■ 

How  defining 

(6.7) 


i(T-l) 


•(T) 


i  a  l(T-l)  h 
\(T)  4  »*h  ” 


(h  /  r  ) 


~*»(T) 


»(T) 


(  /  0  by  choice). 


••(T) 


i(T)  .  1 

*  *  9  ■  n 


(1  i  r) 


i(T)  -  61 
^h  dh 


(h  +  r) 


ws  can  replace  (6.4)  and  (6.5)  by 

(6.8) 


i(T)  k(T)  „  i(T-l) 

he  V  "  h» 


i(T) 


i(T) 


since  a#^  becoaes  coluan  urv  7  •  Cleqrly  only  the  ooluan  x\^  7  and  the  index 
r  •  r(T)  need  be  recorded.  How  Multiplying  both  aaabers  of  (6.6)  on  the  left 
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by  and  on  tha  right  by  we  obtain 


(6.9) 


nJ(T)  h(T-l)  .  J(T) 

%  \  \ 


Equation  (6.9)  ia  tha  haart  of  tha  product  fora  of  inverse.  Applying  it  for 
t  •  lt  2,  . T  and  using  4  to  iodicata  duaey  indices ,  wa  obtain 


(6.10) 


i(T) 

■  “  nK  n> 


V 


4(1) 

*h 


h  Vl  t-2  "t-1 

Hanca  an  aquation  lika  (6.3)  implies  tha  recursive  generation  of  ita  right 
amber  by  uaing  tha  fora  of  given  by  (6.10).  This  ia  easily  dona  as 


follows,  using  (6.3)  for  an  axsapla. 
Let 


fora 


i  -i 

*a(T)  *  aa(l) 

1(1)  ^ 

%  \(1) 


\(2) 


(6.11) 


i(t-l)  -h 


a1 


•(t-1)  *s(t) 


i(T-l)  -h 

-i 

^h 

a (T-1) 

’  “a(T) 

is  easy  to  sea  that 

-i  .  -i 

— r 

i( t— 1) 

*.(t) 

4  °a(t-l) 

V 

(6.12) 

-r 

r(t-l) 

°.(t)  • 

'•(T) 

(i  /  r  -  r(t-l)  ) 
(r  -  r(t-l)  ) 


1(T) 


An  equally  simple  rule  suffices  for  aultiplying  *hv  '  by  a  row  vector  on  tha 
left*  Suppose  it  ia  necessary  to  ooapute 
(6.13)  ft  ^  • 

Va  use  the  transformation  vectors  in  reverse  order  to  that  of  (6.11)  i 
'  Let 


f  - 

xi  *h 

-d)  h(T-l)  .  -<2) 
*h  ^1  *i 


(6.14) 


fora 
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(6.14  cont.) 


Wt)  h(T-t) 
*h  ’'i 


.  7[UV 


-(T-l)  h(l)  .  ;(I) 
*h  ’’l  *1 


Again  it  it  U17  to  •«•  that 

^t.i)  .  ^(t,)  (1  ,  Km)  } 

5<U1>  •  ^  n1^5  (r-r(T-t)) 


(6.15) 


VII  -  THE  PRICIMQ  OPERATION  (Chooalm  Index  a) 

Tha  raw  rector  in  (6,13)  ia  called  tha  pricing  rector.  In  -the  almp- 
laat  problaa  (Phaaa  II,  maximize  x°,  no  infeaalbllitiea),  f^  ■  fcj  •  However, 

in  a  typical  Phaaa  I ,  f^  ■  ,  and  in  ganaral  f^  •  1  for  several  1. 

•4 

Tha  pricing  ractor  ia  appliad  to  a^  , 


(7.1) 


n(T)  a1 
ni  aJ 


(Tha  ara  what  ara  callad  (♦)  ( c ^  -  *j)  in  tha  original  aimplax  method.) 
To  ohooaa  **(1)  Into  tha  baala,  taka 

(7.2)  d<T)  -  aln  4^  <  0  . 

(If  tha  minimum  ia  not  unique,  tha  firat  auch  indax  a  ia  ratainad. ) 

If  all  d^  >  0#  tha  phaaa  ia  ooaiplata.  In  a  Phaaa  I,  thia  la  Terminal  lj 
in  a  Phare  II  it  ia  callad  Taminal  2  and  ia  tha  point  at  which  ona  usually 
axpacta  to  arrlvs  and  quit,  i.a.  tha  optimal  aolutlon  la  attained. 

0 


VIII  -  CHQOSIW3  THE  BASIS  fiOLggt  TO  BJt  REPLACED  (Indax  r) 

Let  tha  currant  baa  la  be  the  currant  optimizing  form  be  row  p 

i(T-l) 

and  tha  ourrant  aolutlon  vector  be  p  '  ' ,  i.a. 


(8.1)  bh  . 

Wa  will  aasuma  tha  aolutlon  ia  faaalbla,  that  la 

(8.2)  0i(T“1)  >0  for  i  >  q 


a 
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where  row  q  it  the  first  actual  restraint  aquation.  Normally  q  -  p*i . 

Tha  usual  orltsrion  for  choosing  r(T)  is  to  chooss  ©p^  as  follows. 
Lst  A  b#  ths  sst  of  indicss  i  >  q  for  which  a1/.,  >  0.  Than 

(3.3) 


V(T) 


■in 
i  £.  A 


'  1 
a 


IV  T) 


>  0  . 


Ths  proble*  of  degeneracy,  i.s.  aultipla  values  of  i ik  for  whioh  -  0,  is 

4is regarded  except  for  ths  following  rule  which  has  proved  effsotive  for 
rsduoing  round-off  error.  (It  ireidentally  breaks  ths  tie  in  Hoffman's  examples 
of  "cycling."  )  If  "  0  and  r(T)  is  aabiguous,  ohooss  r(T)  so  that 

as(T)  largest  possible  (positive)  valufc.  In  case  of  further  ambiguity, 

take  ths  smallest  such  index. 

We  will  modify  (6.3)  slightly.  Lst  R  be  ths  sst  of  indices  i  >  q  for 
which  and  have  ths  same  sign  with  t  0.  Then  lst 


(8.4) 


>(T) 


■in 

ilR 


,  *»<« 


>  0  . 


Note  that  (8.4)  gives  the  same  result  as  (8.3)  as  long  as  (8.2)  holds.  Degen¬ 
eracy  is  resolved  in  ths  same  way,  that  is,  if  0#  take  >  0  and  sax. 

If  no  can  be  chosen,  that  is,  all  ratios  are  non-positive  and  zero 

ratios  have  negative  denominators,  then  has  no  finite  maxima.  A  class  of 
solutions  exist  as  follows i 

(8.5)  mJ(W)  Oh(T*1>  *  e  S(T))  *  9  »i(I)  •  b1 

with  the  value 

(8.6)  pP(T*l)  .  q  *•  • - >  . 

This  is  Terminal  3.  It  cannot  happen  in  a  Phase  I  since  clearly  zero  is  an 


upper  bound  for  the  variable  x 


n*l 
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II  -  5DKMAJBY  OF  CYCLIC  OPERATIOK3  (One  Iteration) 


(9.0) 

(9.1) 

(9.2) 

(9.3) 

(9.4) 

(9.5) 

(9.6) 


(9.8) 

(9.9) 


Test  for  and  of  a  Phase  I  or  for  arbitrary  halt  (external  switch.) 
Datamina  valuaa  of  f^  (dlacuasad  furthar  below.) 

Fom  if)  •  fh  by  (6.14). 

Compute  d^  -  a*  and  chooaa  d^7^  •  min  <  0  or 


taminata  if  all  d^7^  >  0. 


Coaputa  a;/)M  -  a^  by  (6.U). 


i 

‘•(T) 


Chooaa  by  (8.4)  or  tarminata. 


Coaputa  by  (6.6,7),  tranafora  p^7”1)  to  by  ona 


•tap  as  in  (6.11)  and  raoord  s(T)  for  J 


(T) 

r(T)  ‘ 


(9.7)  (Optional)  Print  raaults  of  ltaratlon 


Condansa  and  stora  and  r(T). 

(Optional)  Chack  solution  and  print 


I  -  THE  COMPOSITE  ALOORITHH  (Fomlm 

Suppoaa  that  a  basic  solution  has  baan  obtained  which  is  lnfaasibla, 
that  row  p  is  tha  currant  opt  ladling  fora  and  that  q  is  tha  saallast  row 
index  of  tha  actual  raatralnt  aquations. 

(10.1)  •  b1  ,  h£F  if  h  >  q  and  p^7'1)  <  °* 

Suppoaa  a  vactor  a^Tj  has  somehow  baan  chosan  to  brine  into  tha  basis  and 
that  (8.4)  is  usad  to  determine  r(T).  Aftar  tha  change  of  basis,  tha  naw 


valuaa  of  0* 

.r. 

(10.2) 

ar(T)  .e  > 

S  ®r(T)  * 

(10.3) 

3KT)  .  pl(T-l) 

>(T)  °s(T) 


(i  /  r(T)  ) 


Now  for  i  >  q  and  1  ^  F,  (10.3)  gives  0*^  >  0.  However,  for  i  if,  thara 


ara 


two  cases. 


(10 .4)  If  1£  P  and  a*(T)  <  0,  than  pi(T)  >  p1(>1) 

(10.5)  If  i£P  and  >  0,  than  pi(T'  <  p1^-1) 

Hanca  by  (10.2)  and  (10.3) *  no  infaaalbllltlaa  ara  craatad.  By  (10.2)  and 

(10.4) ,  soae  infaaaibllitlaa  ara  iaprored  or  removed  altogathar.  However, 

by  (10.5)  some  aay  ba  aade  worse.  To  overcoae  this  difficulty,  consider  tha 
v—  i 

function  X  •  >  g  <0  which  if  a  measure  of  tha  infsaeibility  of  a  solu- 

iTP  " 

tion.  We  wish  to  aaxiaise  it  to  taro  and  hanca  wa  aay  replace  tha  aaxljaisa- 
tion  of  pp  by  tha  aaxiaitation  of 

(10.6)  o  -  pp  ♦  X  <  pP 

provided  o  ia  aonotonically  non-decreasing.  Whan  a  reaches  its  smximua,  if 
X  ■  0  than  pp  ia  maxlaua.  If  X  <  0,  than  pP  is  too  great  and  X  aust  ba  in¬ 
creased  without  regard  to  dacraasas  in  0P. 

# 

Let  f^  -  1  if  i£  P  or  i  •  p,  f^  -  0  otherwise .  Than,  as  can  ba  seen 
from  (9. 2,3, 4) 

(10.7)  d<T>  -  tl  «J(T)  <  0  . 

Consequently,  after  change  of  basis,  tha  new  value  of  o  is  * 

(10.8)  o<T>  •  o<T*l>  -  «r(T)  d<T>  >  a'1*1)  . 

If  the  set  P  is  void,  then  (10.8)  is  the  sane  as 

a 

(10.9)  ep(T>  -  Sp(T-1)  -  9r(t)  d<T>  > 

which  is  tha  usual  foraula  for  tha  change  in  the  aaxlaand. 

Now,  however,  if  all  dv.  >  0,  wa  cannot  claim  Terminal  2  without  check- 

(T) 

ing  to  sea  whether  P  is  void.  If  F  is  not  void  and  all  d^  '  >  0,  it  aay  ba 

a 

because  (10.9)  dominates  (10.8).  In  this  case,  wa  aay  scale  down  f  ,  perhaps 

*  Note  that  if  <  0,  than  p1^  <  0  for  i  l  P  and  i  r(T). 

•  ...  .;  *  k .> 


for  hip 
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to  zero,  until  such  time  as  \  •  0.  If  all  d,*^  >  0  with  f  •  0,  then  no 

J  ”  P 

feasible  solution  exists  since  \  <  0  is  it  i  maximum. 

Even  though  \  •  0,  it  is  some  times  desirable  to  set  f  <  1.  Of  course, 
tolerances  must  be  built  into  the  code  in  testing  d,  >  0  since,  if  it  is 
sufficiently  small  in  magnitude,  it  ought  to  be  considers''!  zero.  Varying 
f  has  the  effect  of  varying  this  built-in  tolerance.  Provision  has  also 
been  made  for  setting  f  <  0  which  causes  pP  to  be  minimised  instead  of  maxi¬ 
mized.  This  is  often  convenient  when  experimenting  with  a  model.  The  value 
of  f  is  entered  on  s  binary  card  subject  to  a  switch. 

If  a  <  pp  is  at  a  maximum  and  f^  /  0,  the  machine  will  set  f  •  0  and 
stop  so  that  other  values  may  be  loaded  if  desired.  If  \  •  0  and  f  •  0, 

the  machine  will  set  f  •  1  and  stop. 

P 

The  variable  0P  is  checked  for  monotonic  behavior  each  iteration,  accord¬ 
ing  as  f  >  0  or  f  <0.  This  test  is  suspended  if  4  <  0  or  when  making  arbi- 
P  P 

trary  transformations  when  the  behavior  of  0P  is  unpredictable. 


II  -  PARAMETRIC  PROGRAMMING  (PLP) 

Provision  can  be  made  for  parametrizing  the  right  hand  side  as  a  linear 
function  of  6  >  0,  i.e. 

(11.1)  a*  xJ  -  bl  ♦  0  c1 

where,  if  a  Phase  I  was  used,  c*  must  be  formed  in  the  same  manner  a*  fc  in  (5*26). 
To  do  this,  first  find  an  optimal  solution  for  0  ■  0,  say 


(U.2) 


ui(T-l)  ph(T-l)  .  bi  .  bi(T-l)  # 


Let 

(H.3) 


i(T-D  h 
c 


i 

T  • 


Now  using 


i 

T 


in  place  of  a 


1 

•(T) 


in  (8.4),  a  value 


40  can  be  found  such 


that 


P-842 
P*  25 


with 


and 


(U.M  ,‘(T'1)  (,*«>  .  8r(T)  Th  )  -  b1  .  9  o1 


(U.5)  o1(T*1J  -  #r(T)  r1  >  O  for  1  >  q 


er(T'1)  -  «r(T)  rr  -  o  (r  •  r(T)  ). 


Tht  parameter  6  cannot  ba  ioeraaaad  by  sort  than  with  tha  baaia 

without  violating  (11.5)  but  (11.4)  ia  an  optimal  faaaibla  aolution 
to  (11.1)  for  thia  value  of  d.  Lat 


(11.6) 


b‘(T)  .  „1(T-1)  .  e  el 

r(I) 

1(T)  .  l(T-l)  „  1 

9  S  ®r(T)  1  ‘ 


To  incrM..  0  furth.r,  ^  oust  b.  r.aor.d  froo  th.  b..i.  and  r.pl.c.d 

i  i(T) 

with  aoaa  to  *orB  *  new  80  that 

(11.7)  pi<T>  Sh(T)  -  b1'1) 


ia  alao  an  optimal  faaaibla  aolution  to  (11.1)  for  tha  aaaa  0.  Tha  whola 
procaaa  can  than  ba  rapaatad. 

Tha  indax  a(T)  la  datamlnad  by  tha  rula  uaad  in  tha  dual  aiaplax  algo- 
rithai.  Lat  0  ba  tha  aat  of  indlcaa  j  for  which  a*J  <  0  "  r(T)  )• 
Than  chooaa  by 


lain 


(U-8>  ♦.(!)  *  jit) 


f  ^(T-1)‘J 


-  n. 


j  ; 


>  0  (r  •  r(T)  )  . 


Nota  that 


,(T) 


f®(T) 


-a 


«(T) 


and  (11.8)  ia  tha  analogue  in  tha  dual  problem  of  (8.3). 

If  tha  aat  D  la  void,  than  0  la  at  a  maximum.  If  all  y*  <  0  in  (11.3)# 
than  0  can  ba  incraaaad  without  bound.  Thaae  ara  tha  only  two  automatic 


p-a*  2 

p.  26 


terminations  in  PLP,  The  following  theorem  is  of  interest. 

Theorem:  If  the  choice  of  r(T)  for  (11.4)  is  unique  end  if  D  i»  not  void, 

then  there  exists  e  finite  range  of  0,  $  Tor  which 

ifT-1^  i 

the  solution  obtained  by  replacing  by  ae(T)  *  w^#r*  •(T) 

determined  by  (11.8),  is  both  feasible  and  optimal. 

A  proof  can  be  found  in  Reference 

A  separate  control  code  for  the  oomputer  is  used  in  doing  PLP.  It  always 
starts  from  a  prior  optimal,  feasible  solution.  Due  to  (11.8),  the  iterations 
are  longer  than  the  regular  code. 

Ill  -  MULTIPLE  PHA3£S 

It  is  somewhat  difficult  to  parametrise  the  optimizing  form  since  the 
analogue  of  y*  would  be  a  row  vector  of  n*l  elements.  As  an  alternative, 
provision  is  made  for  multiple  optimizing  fopme  which  can  be  made  to  differ 
by  finite  amounts  in  any  desired  way.  Of  course,  it  is  not  contemplated  that 
two  such  forms  will  be  drastically  different  since  that  is  equivalent  to  two 
different  problems.  In  such  a  case,  it  would  be  better  to  start  the  second 
one  from  the  beginning  or  from  the  end  of  Phase  I. 

It  is  also  possible  to  split  Phase  I  into  multiple  phases.  This  will 
not  be  discussed  further  since  its  application  is  limited  and  it  generalises 
easily  from  the  discussions  given. 

It  will  be  easier  to  describe  the  use  of  multiple  phases  if  a  specific 
example  is  used.  Let  it  be  required  to  optimize  three  forms  and  to  start  the 
problem  with  a  Phase  I.  Then  the  auxiliary  form  must  be 
(12.1)  xJ  .  x"*1  -  0  (J  -  3.  4,  ....  n). 

The  form  to  be  optimized  first  (after  Phase  I)  must  be 
(12. <0  x2  ♦  a2  x^  •  0  (maximize  x2). 
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Similarly,  the  other  two  forme  to  be  optimised  in  turn  must  be 

(12.3)  ♦  aj  •  0  (maximize  x^) 

(12.4)  x°  ♦  tj  ■  0  (maximize  x^)  . 

The  initial  restraint  equations  will  be 

(12.5)  xJ  ♦  xn4i  -  b1  (i  -  4,  5,  ...»  m|  J  -  3,  4,  ...,  n). 

Thus  the  initial  value  of  q  is  specified  as  q  •  4  (total  number  of  phases) 

giving  p  •  q-1  •  3  (number  of  »•  Phase  Il’e"  ).  The  other  two  parameters 

required  are  z  -  1  (number  of  "  Phase  I  *  m '  * )  and  <?  ■  3  (index  of  etsa  row.) 

A'  the  end  of  Phase  I,  p,  q  and  z  will  all  be  reduced  by  1  giving 

p  *2,  q  •  3*  t  •  0  so  that  x *  will  now  be  maximised  disregarding  x*  and  x^ , 

(12.1)  will  now  be  considered  a  restraint  equation  the  same  as  (12.5) »  and 

the  phase  will  be  recognized  as  a  Phase  II.* 

2 

When  x  reaches  a  maximum,  p  will  be  reduced  by  1  to  p  -  1  but  q  will 
remain  at  3  so  that  x^  will  be  maximized  disregarding  x^  and  x^.  Similarly, 
when  x^  reaches  a  maximum,  p  will  be  reduced  to  p  ■  0  with  q  still  remaining 
at  3  so  that  iP  will  be  maximized  disregarding  x^  and  x" ,  In  other  words, 
p  is  reduced  each  phase  but  q  is  reduced  each  phase  only  as  long  as  z  can 
also  be  reduced.  All  three  must  remain  non-negative,  obviously. 


A  Phase  I  is  terminated  when  the  variable  being  maximized  reaches  zero. 

A  Phase  II  terminates  when  all  d^  >  0.  These  criteria  are  not  always  equi¬ 
valent  even  in  Phase  I,  If  b*  is  representable  with  fewer  than  a  of  the 
columns  a*  (i,  J  >  0),  then  several  artificial  columns  may  remain  in  the 
basis  at  zero  level  at  the  end  of  Phase  I.  In  this  case,  the  d^  corresponding 
to  the  Phase  I  pricing  vector  will  not  all  be  non-negative.  If  Phase  II 
starts  with  artificial  columns  in  the  basis  (other  than  a*  and  a*^),  then 
a*^  may  be  eliminated  in  Phase  II  but  one  a^  will  always  remain. 


•J 


?•$& 

n.r* 


Rrfoien^cs, 

The  " oi  lman-  rro'”'  >  material  on  the  ’*ev:‘ aed  simplex  method. 


1.  D'intx'r,  t.  T>.,  Ale::  f  ••der  inf  Pill  >  "olfr ,  'Ves  on  Ilr.o-'r  rropvaamifte.  Part  I 
?hr  raxy  i.rir  r.iMir  ?g7nor>  -or  r:i::T?7'  irn  a  ijitej  very.  ui;de--  liheap. 

: •%  ui it.  t.-t  •■—"•.  x-r  n,  tv  east  ^orr.,  r.antn  KonK,  raii-'. 
_ , _ ,  r  v  n,  r>u\ij?y  theo;  gr, ,  rm-i?6$,  The  pajt  "crp. 


",  Dan ta 0.  ?),,  V.’m.  '  r'har  ’-K:*  n ,  !’o* vc  on  Linear  Pro,~ran»nlnr«  Part  V.  .■.LTST^u.^ 
,.ir.oer?ig.-  :hk  V3E  ysvrrn-  TtEfM  >r>  ut  :nr  a  Produr-  ror»  for  the 

Inverse  ,  ■■' t  Tlic  J.l?  Oo'-p. 

.-  n,  •*».,  !  ..r-:r,:-rtri),  h2vtiu.?ie:v  uai  py.TET^iors  qr  the  rmsn  pimply 
tariia-.  1  ■ :  ,  ^  ~<r  . 


The  foil  curin'-  d'r-M.rr  vjy'iv.l  rinplor  method. 
r.  ro».loe  ronn’er  on  Kor.n'-’-nh  *>-.  !*,  AffTrVTTr  AiAIYSi::  CI~  J^OPHTTIO::  AIID  AUPCA- 
TICTI,  ?.  J.  "ooTKiunr, ,  -  ’-‘‘or.  h-n.  Ycr-h,  John  .’lie;-  and  2ow,  ino.,  1X1» 
(fee  espe-  'all:  'Ttaptere  JEI,  DI,  miT) 

6.  "V.mec,  A.  one  17,  W.  "ooper  and  A.  Henderson,  AT?  ITEBRODIFTIOK  to  LH3SAT  l?p- 

CHAM-ITTIS ,  ITev  York,  JoV  <’iley  and  Tons,  In-.,  l?r>5. 

7.  Eiromann,  ,  I.IKR\i:  PR(»TRA!lMIITn ,  Oaarherly  of  Applied  Mathenc'lrs,  Vol.  XT  IT , 

I!o,  Oc  .ober  l?r<% 
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CJ'H ‘j.'I*’.' 7TT’  Cv  Tin?  ''•^TT 


The  date  assembly  -ode,  the  none  ’m]  <■?,  1.*  ar  .isuembiy  ^rorr-n  ami 
contains  the  ubup.1  ~anvorBion  routines,  etc.  as  veil  as  several  roeeini  eub- 
roet.'nee  denlmed  fer  4 he  part*  *ular  ■'orvertlonn  used  Jr  ‘hie  system.  A  de- 
tabled  description  o*'  its  orran* nation  voOl1  serve  no  -  refill  purpose  here. 

Complete  1  in4  >rs  of  the  rrormn  arc  available  for  ‘hose  -.-ho  have  a  technical 
‘nt erect  in  the  codinr.  The  preot.ra+ ior.  o  '  the  inputs  are  described  under 
IIIRfr  and  the*  output  is  C's-ussod  *  nder  CPF!  /TIOT:  OF  THE  COPER.  We  v’ll  mere¬ 
ly  note  here  *hat  data  is  processed  either  f-am  -ar >■  or  from  tape  '  and  that 
output  io  to  tape  ‘  and  prnchi  .  -ar>ir.  (tlnar:.)  Jin-e  the  entire  pro rrori  <c 
too  lonr  for  the  available  part  of  core  ctor.are,  a  part  of  the  data  assembly 
pro -ran  is  stored  on  drum  1  (automat  ally,  dtirinr  load inr)  and  recalled  to 
core  Sworaye  vhor.  n-.-'-ad . 

The  'or'i — r  -r*'  ir  v’>'-  o  r  oryoaisa'^on  of  the. mu* n 

routiner,  ur  :os^  •  t  e  >.  c*  . 

Coneoptusll/,  tho  hlgh-opM d  iton  (IBS)  1*  dlvldod  into  two  Min 
Motions t  (1)  tho  progroao,  eo no touts,  poroaotoro  and  temporary  otorogo, 

(11)  tho  dots,  both  original  had  transformed. 

SlalUrlj,  auxiliary  otorogo  Is  dlvidod  Into  two  parts,  ono  port  os  poimon- 
ont  otorogo  for  (1),  tho  othor  so  permanent  otorogo  for  (11).  Thus  the 
ontlro  HSS  Is  la  o  osnso  ^temporary  storogs"  .  In  proetlco,  o  consldoroblo 
port  of  ooetlon  (1)  of  HSS  rsMlno  In  toot  throughout  o  run  but  oaa  bo  ro- 
otorod  from  auxiliary  storsgs  (drums)  If  necessary.  Ths  odvant  of  tho 
extremely  rolloblo  coro-otorsgs,  however,  makes  tho  nsod  of  rootorlng  unlikely. 

Vhon  HSS  roochoo  tho  sisss  now  contemplated  (that  lo,  32,766  words  of  ears  otorogo), 
suxilisry  otorogo  for  ths  cods  and  for  cortoin  dots  will  bo  onooeoosory.  Tho 
aoln  uso  for  auxiliary  otorogo  lo  for  tho  a*  natrlx  and  tho  vac  tors. 

Hognstle  topos  oro  uosd  for  both.  On  tho  IBM  704,  throo  topos  oro  set us 11/ 
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um4  for  the 


rectors 


as  explained  below. 


The  data  section  of  K3S  la  dirlded  into  four  regional 
the  H-region  for  the  baeia  headings:  m*l  words 

the  V-region  for  the  values  p*  of  the  current  solution,  aain- 
tained  in  double-precision,  floating  point  i  2(»*1)  words 
the  W-reglon  for  work  spaoe  in  generating  a*.  ^ .  and  and  other 

purposes:  2(a*l)  words. 

the  T-reglon  for  tsaiporarily  holding  the  a*  Beatrix  or  vectors 

or  as  such  as  possible  of  either:  the  remainder  of  section 
(ii)  which  should  be  at  least  4(*+l)  words.* 

The  sis#  of  these  regions  is  s  function  of  the  number  of  restraints.  Their 
origins  are  computed  by  the  data  assembly  program.  in  PL P,  ?- region  is 
used  as  a  second  W-region  in  pricing.  A  duplicate  of  K-  and  y-reglone,  as 
wall  as  the  original  b*,  is  kept  in  auxiliary  storage.  This  allows  re-starting 
after  an  error  and  checkirj  a  solution  tqr  computing  and  printing 


.  ,!<»)  ,»<’>  -  b1. 


For  PL P,  it  is  also  necessary  to  keep  c*  and  T1  in  auxiliary  storage. 

The  prograa  section  of  HS3  is  divided  into  s  even  regional 

(a)  Temporary  storage  called  COMOV. 

(b)  The  aain  control  region,  oalled  the  NCR. 

(c)  A  sub-routine  for  doing  double  precision,  floating  point  addition, 
called  DPFADD. 


(d)  A  similar  sub-rout  ins  for  multiplication,  called  DPFWUL. 

(s)  A  routine  called  the  distributor  (DISTRB),  explained  below, 

(f)  Spaoe  for  the  largest  sub-routine.  The  first  location  is  called 
SRORIG, 

*  The  origin  of  T-region  has  been  arbitrarily  set  at  2048  for  all  problems. 


(g)  Unlramal  coast  ant a,  pinwUn,  origins  ate.,  oallad  K-rtglon. 
Tha  dn4  for  OOMfOl,  DPFADD,  DPFNUL  and  K-ragion  la  obrioua.  Tbolr  looa- 
tlaas  art  ptraaasntly  f lxtd  and  may  bt  rtftrrtd  to  fret  any  program  In  tho 
systsa.  Ail  sub-routlntt  art  olottd.  Tht  K-ragion  eontalat  otrtaln  calls ' 
whota  coot ant a  art  fixsd  only  for  part  of  an  1  to  rat  ion  or  othtr  sub-ssqusncst 
of  tbt  pro  bits.  Thtrt  art  conrantiona  oo  whan  and  by  what  routint  tho  at  art 
to  bt  ebangtd.  OOMNOV  it  always  availablt  to  any  routint. 

Tht  function  of  tbt  NCR  it  to  aaka  tht  aajor  doc  1a loot  and  call  for 
tht  proptr  stqutnot  of  optrationa.  Host  of  tht  actual  work  it  daltgatad 

to  aub-routinas .  Bttidtt  tht  NCR  for  tbt  at  in  algorithm,  thtrt  it  ont  for 

> 

PL?  and  ont  for  ra-lnwarting  a  batia.  Othtr  auxiliary  programs  art  tatily 
ertattd  by  coding  a  now  NCR. 

Tbo  NCR  oalls  for  a  aajor  opt  rati  on  by  linking  to  DISTRB  with  a  pssudo- 
oporation.  DISTRB  calla  in  tho  proptr  tub- rout  lno  frea  auxiliary  t  tor  ago 
and  loadt  it  into  KS3  atarting  at  SRORZO.  Control  it  than  tranofarrod  to 
tht  sub-routina  which  roturnt  control  to  tht  NCR  after  completing  its  func¬ 
tion.  If  soma  othtr  trrangtasnt  for  linking  to  tub-rout  into  ia  das  1  rad  (at 
for  tBGMpla,  whan  HSS  is  rory  larga)  it  ia  only  nacssaary  to  ohango  DISTRB 
to  arrangs  for  tht  proptr  routint  to  taka  osar  in  tht  proptr  way. 

Thtrt  art  fourttan  standard  aub-routinaa.  Othart  may  bt  addtd  for 
tptclal  purpoatt  if  da  el  red,  up  to  tht  limit  of  auxiliary  atoraga  to  hold 
than.  Nora  importantly,  if  an  lap  ro  rad  or  a  aptclal  rtrslon  of  ona  it 
dawalopad,  it  can  rtplaot  tha  old  ont  aartly  by  axe  hanging  tht  proptr  binary 
cards.  Thttt  aub-routinaa,  togathar  with  DPFADD,  DPFNUL,  K-raglon  and  a 
aptclal  loading  routint,  axlst  in  binary  cards  which  conatltuta  t  basic  dack. 
Tht  origins  card  produetd  by  tha  data  at  stably  la  put  ahaad  of,  and  tha  NCR 


behind  this  basic  deck. 


As  soon  as  the  special  loading  routine  is  in  KSS,  it  takes  control 
and  loads  the  remainder  of  tbs  basic  deck  and  the  NCR,  storing  than  in 
auxiliary  storage  (druas).  Since  many  of  the  sub- routines  require  addressee 
to  be  set  which  depend  upon  m,  a  loading  interlude  is  performed  to  do  this 
initialisation.  This  is  accomplished  by  an  initialising  routine  which  goes 
with  a  sub-routine,  the  two  being  loaded  simultaneously  into  K33.  Control 
is  then  sent  to  the  initialising  routine  which  does  its  job  onoe  and  fbr  all 
*  and  returns  control  to  the  loader  .  The  loader  stcree  the  initialised  sub¬ 
routine  in  auxiliary  storage  and  builds  a  catalogue  of  locations  into  D1STRB. 
The  NCR  takes  over  control  when  loading  is  complete. 

The  fourteen  sub-routines  (called  «  codes*)  are  as  follows.  Some  are 
used  for  multiple  purposes  which  are  controlled  by  internal  switches. 

Code  1.  Load  the  binary  oards  produced  by  the  data  aeeembly  and 

•tore  in  appropriate  placet.  If  the  a*  card*  are  included,  they 

J  • 

are  transferred  to  tape.  Subject  to  a  switch,  a  value  may  be 
loaded  from  a  special  card  for  f^.  Ry  meant  of  an  internal  switch 
and  the  use  of  control  cerds  which  it  loads,  this  load  routine  is 
able  to  distinguish  various  casts,  initial  start,  re-start,  start 
of  PLP,  re-start  of  PLP,  so  that  the  proper  information  oan  be 
stored  in  auxiliary^ 

Code  2.  Store  the  current  solution  in  auxiliary  (drums).  The 
"  current  solution"  is  defined  as  K-,  H-,  and  T-reglons. 

Code  3.  (Normal)  Form  the  row  f^  in  W-region,  and  record  whether  any 
variables  are  infeasible. 

(During  PLP)  Fom  1  in  V-reglon  and  ‘I  ln  W-region. 
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Cod* 


4.  (Normal)  Compute 
(During  PL P)  Compute 


«hT)  •  fl  'i(T'1>  l»  w-r«glof\. 

in  V-region  and  n^”^^in  W-r*gion. 


Cod*  5*  (Normal)  Price  the  matrix  aj  and  choo**  dg^. 

(During  PL P)  Cho***  11  described  in  Part  A ,  Sect.  XI. 

Subject  to  an  external  switch,  recognize  "  Curtain"  marks  in  a*  . 

Subject  to  instructions  from  the  MCR,  make  the  following  check t 

(Normal)  If  d<7)  +  0,  (during  PLP)  if  n^7”1'  aj^  0  for  J  ^  jj.Jjj* 

then  test  for  J  occurring  in  jjj7”^ .  If  it  do**,  an  error  has 
ooourred.  This  check  is  very  valuable  for  detecting  errors  before 
much  more  computing  is  done.  It  does  consume  some  amount  of  time, 
however,  which  increases  with  n. 

Code  6.  Load  from  a*  matrix  into  specified  region  as  double- 

precision,  floating  point  veotecr. 


Code  7.  Multiply  nj^7"^  a*(T)  ^0r  tJn^  ▼•otor  in  specified  region). 


Code  8.  (Normal)  Choose  the  index  r(T)  *s  discussed  in  Paz*t  A, 
Sect.  VIII.  For  arbitrary  transformations  or  inversion,  choose 
r(T)  by: 


r 

MX 

i 

*»(T) 

"  ia 

“•(T) 

where  A  is  the  set  of  indicee  i  for  which  (artificial). 

Code  9.  (Normal)  Compute  from  and  transform  0^7-i^  to 


0 


KT) 


(PLP,  first  entry)  Compute  0 


i(T) 


and 


©r(T)  from  0i(T‘1),  r(T)  and 


Yi(T“1).  (Part  A,  S*ot.  II).  Then  compute  b1^  -  b1^7*1^  e^c1. 


#  “part  A"  is  used  here  to  refer  to  DETA IIXD  JR  1TE-UP . 


l(r)  i  • 

(PLP,  second  entry)  Compute  r\r  '  from  and  transform 

Y^^~^  to  (All  thess  operations  are  very  similar,  even 

■ore  than  appears  at  first  glance.  The  variations  are  Merely 
•witches.) 

1(T) 

Code  10.  Delete  seroe  fro*  r  v  '  and  index  non-sero  eleaents.  Store 

r 

condensed  vector  in  auxiliary  storage  (druias;  cf.  Code  13). 

Code  11.  Multiply  out  -  b*  -  £*  taking  fro* 

original  a*  Matrix  by  referring  to  jj^.  £  ^  is  left  in  T-reglon 
and  b*  in  W-reglon  for  printing. 

Code  12.  Print  prograM.  This  program  is  it  a  elaborate  and  longer 
than  the  other  codes.  Special  provisions  for  it  need  not  be  dis¬ 
cussed  here  incept  to  say  that  appropriate  captions  and  identifi¬ 
cation  of  columns  are  printed  for  each  type  of  print-out,  of  which 
there  are  14.  The  print  output  can  also  be  put  on  the  fifth  tape, 

if  desired,  for  later  printing.  In  all  print-outs,  there  are  5 

(T}  i(T) 

columns  of  which  the  first  three  arej  »  0  '  and  1*0,  ..., 

The  last  two  are  the  contents  of  W-  and  T-regions. 

Oode  13.  This  performs  an  »%nd-of-stage "  it  which  soae  of  the 
are  transferred  from  drum  to  tape.  Gee  below. 

Code  14.-16.  (Undefined  ,  except  for  code  1^  used  vith  a  special  . ) 
Code  19.  Automatic  restart  program  for  recovering  after  an  error  by 
returning  to  the  beginning  of  the  iteration.  Clearly,  this  is 
highly  dt pendent  on  the  psrtlcular  machine.  The  important  points 
to  note  are  that  ail  programs,  the  original  data,  the  current 
eolutlon  at  beginning  of  the  iteration,  and  the  r,  '  1  rectors  must 
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bs  intact  aomewhsre  tn  the  raaohina.  Also  trds  code  must  have 
been  previously  instructed  where  to  return  control  to  and  he  able 
to  restore  all  equipment  (e.g.  tapes)  to  the  proper  poeltlor.e. 

It  suet,  of  course,  be  activated  by  tom*  external  means. 

Besides  the  physical  organization  above,  Ujere  are  also  dynamic  groupings, 
several  of  which  have  already  been  discussed:  the  iteration,  the  phase, 
maximisation  of  A<0,  arbitrary  trans  format  ions,  etc.  Another  grouping  is 
called  s  stats  and  alwaya  consists  of  an  integral  number  of  iterations.  It 
has  nothing  to  do  with  phases  or  other  mathematical  aspects  but  is  simply  an 
operational  device.  It  wss  devised  for  the  IBM  701  and  has  been  carried 
over  to  the  704,  for  ellghtly  different  reasons.  Some  analagous  arrangement 
la  probably  necessary  on  any  machine. 

Most  of  the  advantage  of  condensing  the  vectors  would  be  loat  if 

a  separata  acosas  to  auxiliary  storage  had  to  be  made  for  eacn  one.  It  la 
desirable  that  aa  big  a  chunk  of  ther«*  transformation  vectors  as  possible  be 
recalled  at  one  time.  The  limiting  factor  is  the  size  of  T-re^ion.  Hence 
aa  tha  are  generated,  they  are  stored  on  a  drum  until  or.e  mor'  would 

exceed  the  capacity  of  T-reglon.  At  this  point  an  end-o f-stage  procedure  ie 
performed.  Before  describing  this,  it  is  necessary  to  explain  the  use  of 
three  tapes  for  the  vectors. 

The  tapes  on  the  704  can  be  baok-apaced  and  additional  records  can  be 
added  to  an  axisting  fila.  However,  they  cannot  be  read  in  a  backward 

direction  and  the  back-space  is  somewhat  slow.  Hence  the  records  are  stored 

(T) 

in  reverse  order  on  a  second  tape  for  use  in  computing  .  The  third  tape 
is  used  alternately  with  the  second  from  stage  to  stage.  Though  thi->  costs 
somt  copying  time,  it  doss  provlds  a  apart  (and  chscked)  copy  of  the  tApe  at 


i 


I 


o; 


The  end-of-etage  procedure  Is  u  fo^.i  ws  •. 


(1) 

(2) 

(3) 


Compute 

Print  j[T),  0i(T),  i,  b1,*  *. 

Transfer  the  on  dnsn  to  an  aciitional  record  on  tha 

forward  tapa. 


(4)  Transfar  tha  on  drum  to  tha  first  racord  on  the  fraa 

backward  tapa. 

(5)  Copy  old  backward  tapa  to  raaalndar  of  new  backward  tapa. 

<  t  '  *  r  t  ^ 

(6)  Punch  out  binary  cards  containing  K-region,  J’  ;  and  p  '  '. 

During  PLP,  also  punch  b1(T>  rl(TK 

(7)  Adjust  tha  nacassary  bookkeeping  parameters  u 

An  end-of-etage  aay  ba  forcad  by  an  external  switch.  It  also  occurs  s*  tha 
and  of  Phasa  I  and  terminatione .  nteps  (1)  and  ( j  above  aay  be  forcad 
without  tha  others  at  tha  end  of  a  cycle  by  \n  external  switch.  This  would 
ba  desirable  aoove  tape  manipulations . 

To  restart  from  the  and  of  a  stage,  it  is  only  necessary  to  use  tha 
punch-outs  to  replace  the  corresponding  original  data  cards  and  put  tha 
tape*  back  on  tha  same  units. 

Similar  simple  hand  collating  of  blocks  of  punch-outs  with  the  original 
data  cards  (binary)  and  changing  of  ths  MCR  ar#  all  that  are  necessary  to 
aet  up  the  deck  for  PLP  or  re-inverting  the  basis.  Special  short  MCfi's  are 
usually  accumulated  for  such  things  as:  transforming  a  new  right  hand  side, 
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,o*ratla«  40d  pristine  tb.  dJf  »lt.rln«  «  Ttctor  In  *h.  b*«U,  *tc.  * 

reutlna  1.  pwlM  f«  cpyln.  .  ^  Up.  in  r.«r«  oni.r  .o 
ihfti  a  prohlM  can  to  picked  up  after  tod  luck  with  tape*. 


M4a 

1.30 


The  data  mmII j  prsgran  (which  will  be  discussed  only  briefly) 


as  identification  card  followed  by  a  card  with  tha  5  peraaeteri.  Aftar  this 


tba  basis  headings*  If  say,  FIRS?  B,  tha  b  ,  ato.  aa  daacrlbad 


ba  ara  sntsrad,  with  thalr 


DM,  Only  aoa-ssrs  elmmrts  ara  entered,  with  thalr  proper  ladloM.(Orp 
will  canes  bo  difficulty,  aaly  a  waste  of  tlas  aid  apaaa.)  Tba  yr#» 


pwnbae  la  blaary  oarda ; 


Tha  ldentlfloatlan  card  (blaary  cadad)  which  amt  ba  placad  after 


tha  mttrlx  traacltlcB  card 


lead  lac  by  sals  oede. 


A  aalf  lead  lac  eorlglM  card  ndeded  by  tha  lead* 


tha  aalm 


tinea  which  anat  ba  placed  a*  tha  front  of  tha  mim  oada  deck, 


fWa  blank  separators,  (discard) 

Initial  X*,  B-,  and  T-reglane  fallawad  by  a  blank  (retain.) 


Right  hand  aide  for  lead  lag  In  V -region. 


(Optional)  Blade,  auxiliary  right  hand  aide  far  loading  In  V«reglga» 
(optimal)  Blank,  natrlx  far  landing  la  T-reglcn.  Multiple  rscnriu  me 
separated  by  blades  (retain  but  net  aftar  last  reoerd.) 


Hgtrlx  transition  card  (9  punch  In  col. 57  only.) 


Ter  srrsr  detect  ion,  eoa  list  of  stops  - 


Tha  erderlng  of  tha  cards  la 


ldpertant.  XT  ♦< 


oarda  ara  used,  they  Boat  net  ba  planed  after  the  last  card  af  a 


they  lead  late  the  highest 


call  used  by  tha  reatlas.  She 


ft 

A  *J 
r  . 


.V 

>: 


q 


4m 


l  -  na 


track  af  first  and  last  lecatlm  af  each  /eotiae  except  far  MB'e,  and  uses  tk>  11  J 


itlm  for  loading  tba  rautlnaa  on  tba 


and  building  tba  table  for  non 


«  x  v 


by  yiPTU.  Following  each  routine,  a  title  eari  Idagtlfiea  the  cede  Jmt  lddd«‘!3j^j 


The  9  row  left  la  blade.  Tba  9  m*  rlgrt  contains  tha 


im  t»  <*•**.  -4 
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■nit,  •db-JCl'a  are  Identified  by  taro  and  MSB's  by  minus  taro  In  9  rov  right. 
The  print  jrograa  la  loaded  on  drtm  2  first  as  though  it  vara  a  aub-ITK.  Bov- 
war,  this  oausss  no  difficulty  alnca  DISH®  doaa  not  racall  the  print  program 
to  BB8  directly.  A  pseudo-code  12  (loaded  in  the  normal  vay)  recalls  the  real 
print  program  frcn  tha  knovn  location  at  the  beginning  of  drtm  2.  The  print 
program  occupies  tha  high  part  ef  T-regioo  vhan  in  use.  If  another  rub-ICH  la 
loaded  later  in  the  loading  sequence,  DI8TKB  recalls  this  vhan  Coda  0  la  called 
for  by  the  MSB.  DXBBLB  dees  not  reeagnlsetode  -0”. 

The  binary  cards  a re  id  amt  If  lad  by  Boiler  1th  character*  In  cola.  73-00. 

The  basic  deck  la  ae  felloes. 

CIS  KBS  1  card  self -loading  program  to  claar  Index  rag  later*  and 

load  nart  card.  (Preeadad  by  erlglne  card  punched  by  data 
aaeert> ly  vhich  laavaa  Index  rags.  in. roper ly  far  BIA&81XR.) 

PAB8IXB  1  card  aalf-loadlng  lead  program  far  loading  LPBABX. 

UTASBOl  Several  service  rout  Inset  EEBOOT,  a  drum  n  bootstrap"  for 

.*  activating  Coda  19;  DFTA JX>j  DTTMJL;  DISH®  vlth  blank  table; 

XJBA8P3  and  SUCID,  the  actual  lead  program  dlacuaaad  above. 

XJBAJBB  Transition  card  to  give  control  to  KHCLD. 

LHIBT01— XJRISTB5  The  actual  print  program, 

lift  URL  Title  card  for  print  program. 

LPCD0101— LPCD0105  Coda  01 

UCB01TL  Title  card 

UCD0201— LPCDOBOC,  UCD02TL  Code  02 f  etc. 

1 JCD0301— ZJPCD0903 »  XJCDOfffc 

UCD0401— : ipcdo4io,  ucDobn 

MCD0501— UCD0511,  LPCD05TL 
UCD0601— LFCD0603 ,  LPCD06TL 
IPCDC701— UCD0706,  LPCDOTO 
1PGDO0O1-- LPCDO0O8,  LPCDO0TL 
1PCD0901— LPCD0909 ,  LPCD09TL 
LPCD1001— LPCD1004,  LPCD1OTL 


LPCD1101— 1PCD1109  $  UCD11TL 


LPCD1200  Pseudo-code  12  for  c  all  U|  print  pragma  Into  BBS. 

IPCD12TL  Title  card. 


UCD1501— IFCD1309,  LKD13TI 
UPCD1901— IPCD190*,  LPCD19TL 


LPBDTT01 

0 

Lwimno 

XJBIWTL 


Bub -ICR  for  inversion  code.  Nay  be  left  in  baeic  deck  if 
desired,  as  long  ae  no  other  cub -ICS  la  ueed  in  which  cate 
the  capacity  of  the  dr\aei  nlgrt  be  exceeded . 


THE  THRU  mm  ICR's 

§ 

The  baeic  ICR  ie  for  the  composite  slag)  lex  algor  it  ha  Including  Phaee  I, 
■ultlple  phaeee  and  arbitrary  transformations.  The  carde  are  identified  by 
LPC0MC1 — LIC0MK16,  IKOMffL. 

Theee  follow  the  baeic  deck  described  shore. 

A  special  ICS  le  ueed  for  PIP.  The  oarde  are  identified  by 


LPPLWDl— LPTUIQO,  LPPLPIffL. 

When  starting  a  ran  to  do  PIP,  these  cards  are  us<»d  instead  of  LFCCMOl,  ate. 

Another  ICR  is  ueed  to  re- invert  a  given  basis.  Put  LPBDJYOl,  etc.  at 
the  end  ef  the  baeic  deck  (see  abort)  and  use  the  following  ICR  carde. 


IPIWTMai— IPTR71Q5,  IPIRTHTL. 

The  sub-ICR  operates  first  to  organise  the  vork.  However,  the  operation  of  this 
is  autoaatlc .  LP0IRTO1— LPBDfnO,  LPBUTVTL,  LPIRTM01— IPITfTK15 ,  IPOfNEL  can 
all  be  left  together  and  considered  as  one  ICR,  if  desired. 


SWTCHPP  ICR’s  DQRXH3  A  ROT 

It  ie  sonet loss  desirable  to  start  with  one  NIP  and  then  switch  orer  to 
another  one  without  getting  off  the  nachine.  It  is  not  necessary  to  use  the 
punch -oute  and  reload  the  basic  deck  to  do  this.  Severer,  if  a  switch-over 
le  to  be  nede  to  the  inversion  ICR,  the  sub-ICR  oust  have  been  loaded  vith  the 
basic  deck  Initially. 

The  special  switch-over  procedure  uses  Code  19,  sense  switch  6  and  a  special 


transit lee  card. 
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There  are  tvo  of  these,  one  for  uie  vlth  UC0M4  or  LPTWM 
■arked  PHOT  TP.  The  other  le  for  uoe  vlth  LFPLFM  explained  below. 

To  evltch  over  to  either  LPC0##4  or  LP1ITVM,  proceed  as  follore: 

1.  Put  the  foil ov In*  sequence  of  cnrdo  In  the  render: 

HCR  without  title  card 

PDTCP  TH 
Blank  cart] 

Matrix  transition  card  (9  In  col.  37) 

Idezrtlf  lcatlon  card  punched  by  data  assembly 
Three  blank  card* 

2.  Put  Switch  6  down. 

3.  Push  CXSAB  button. 

4.  Push  LOAD  BRIM  button. 

3.  Her  Is  read  In  and  Machine  halts.  Put  Switch  6  up. 

6.  Push  START  button  to  proceed. 

If  It  desired  to  read  a  scale  factor  card  (f  ).  then  Switch  6  la  left  down  In  3 

P 

and  the  fp  card  Is  put  after  one  of  the  blanks  following  the  Identification  card. 

To  switch  over  to  LJTLFM,  the  procedure  le  the  sane  except  for  the  cards 
used  which  are: 

LPPIPMD1— 1PPIP10 

The  e1  cords  punched  by  the  data  assembly  (which  load  to  V -region) 
PLP  mCL  FXXSPTlf  (card  is  punched  this  way)  and  tvo  blank  cards . 

An  f  card  cannot  be  loaded  for  PLP  in  this  Manner. 

P 

upoih  or  rm 

The  codes  have  been  designed  to  operate  continuously  until  either 

(a)  a  phase  is  complete, 

(b)  switch  .1  is  put  down  (arbitrary  end-of -stage  and  halt),  (also  sv.6,  see 

bo  lov  ) 

(c)  the  value  of  f  Mist  be  changed,  car 

(d)  an  error  is  detected. 

Whenever  a  stop  occurs,  hitting  the  START  button  will  cause  the  program  to  do 

the  "right  thing"  ,  insofar  ai  possible.  This  nay  be  either 

(a)  to  continue  the  calculations, 

(b)  to  try  again  the  sequence  of  operations  in  which  an  error  was  detected,  or 

(c)  to  lock-up  emd  do  nothing. 


M4  2 

vM 


If  the  chine  he*  performed  everyth  in*  required,  It  vill  leek -up  sith  ell  feur 

im«  llghte  an.  If  the  error  1  s  euch  that  it  cannot  be  correctod,  the  automa¬ 
tic  drvmi  reetart  (Code  19)  nay  be  ueed  to  repeat  the  iteration.  This  is  acti¬ 
vated  by  clear  lag  8ft  (CLSAP  button)  and  pushing  the  LOAD  EFTH  button.  Be  sure 
switch  6  Is  up  shen  doing  this. 

All  stops  are  described  in  the  list  given  later  belcnr.  The  cnee  normally 
to  be  expected  ore  0,  1,  6,  7,  and  27.  After  certain  stops  (e.g.  5,  122)  Svltch 
1  nust  be  put  dovn  before  hitting  STABT  (to  get  an  etkl -of -stage)  in  order  to 
stop  the  nachine  at  the  proper  point.  The  job  nay  be  taken  off  the  machine 
after  any  md-of-stage  and  restarted  from  this  point.  The  punch-outs  replace 
the  correepend lng  original  binary  data  cards,  l.e.  K-,  H-,  T-reglsne  and,  for 
PIP,  V -region  (}  ^).  The  binary  data  cards  for  COMO  and  PXPMCB  are  cosqmrod 

belov. 


cgnggo  or  qg  b  jvm  data  cards 

These  cards  folios  lmned lately  after  the  K7R  title  card. 

tfcox  imsm 

K- region 
H-reglon 
▼-region 


W-region  (b^T*) 

/  T-regicn  (y1^)  * 


Blank  card 

W-region  (c1) 

Blank  card 
i 


Blank  card 
W-repion  (b*) 

Blank  card 

aj  natrlx  (T-regicn),  (Optional) 

Matrix  transition  card 
Identification  oard 
Blank  card 


t  card 
P 

2  blank  cards 


(Optional  Y 


{ 


trix  (T-reglon) 


insert 
to 
start 

np 


Matrix  transition  card 

Identification  card 

Blank  card 

f  card 
P 

2  blank  cards 


P-042 

P>5 


(Key**)  fro*  la  at  page) 

Current  7*^  Teeter  does  not  exist  initially.  Par  starting  PLP  it  ie 
substituted  for  by  a  card  vhlcb  la  blank  except  for  9-punch  In  some  (any) 
celuam  58-72.  Thia  doea  not  replace  the  blank  following. 

**  Read  only  if  Switch  6  down.  The  two  word  a  in  9-row  are  the  value,  the 
binary  point  being  between  column  26  and  27.  If  fp  ia  to  be  negative,  put 
a  9-punch  In  both  colunma  1  and  57 . 


To  at  art  an  inversion  run,  use  aane  data  aet-up  an  for  XPCOMM  with  the 
punch -out  n  vhoee  H -region  npecify  the  basis  to  be  inverted,  everything  is 
automatic  provided  the  dee  Ignat  Iona  tJPOOl,  UP002,  etc.  have  been  used  for 
legit inate  unit  vectora  5^  ,  ,  etc.  If  any  J^<  0  are  left  o  x  “ron 

original  aet-up,  they  will  be  treated  an  zero  and  discarded. 


•  *1  •* 


The  captions  on  the  print -outs  describe  the  reason  for  printing.  The 


identification  card  la  always  printed  first,  followed  by  a  line  with  iteration 


ntntxar  (T),  stage  mok>er  less  one  (l.e.  maker  stages  coagxleted),  form  number 
of  w^****"*  (l.e.  current  value  of  p),  caption  and,  for  sons  prints,  *3  IS  1  * 
with  appropriate  vector  nans  filled  in.  This  naans  a^  just  cans  into  the 
baala  or  is  to  be  multiplied  by  B  far  unbounded  solutions.  The  third  line 
ie  a  set  of  headings  for  the  five  colunna  of  printing  of  which  the  first  three 
are  always  1 

J  for  ;  BfEA  far  ;  I  for  running  index  1. 


The  fourth  coluxn  nay  be  headed  by  any  af  the  follwlng: 
B  far  b1 

I  for 


0 

A 


(or  b1^  In  PIP) 


for  y 


i(T-l) 


for  0 


f(T) 
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p.tt 


The  fifth  colimm  my  bt  either 

M  tar  £ *  (error) 

,(*) 


PI  for  * 


or,  In  one  Instance,  blank. 


or 


The  fourteen  caption*  vlth  the  correspond  Inc  head  Inga  for  the  fourth  and 
fifth  coltamia  are  as  follows: 


1. 

CTCLX  PRIFT 

S  IS  _ 

1 

PI 

(forced  by  Switch  k) 

2. 

HIV  SOUTH  OF 

S  IS  _ 

B 

PI 

(occurs  every  Iteration  in 

FI#  unless  $r  -  0.) 

з. 

и. 

CHICK  SOUTH  OF 

TO  OF  FHA8X  OKI 

B 

B 

Et 

W 

(forced  by  Switch  6) 

5. 

TO  OF  STAG I 

B 

n 

(dnm  3  full  or  Switch  1  down) 

6. 

10  FIASIBLI  80U7H0K 

B 

n 

(see  stops  2  and  b  below) 

7. 

FZASIBLI  80UTH0F 

B 

m 

(see  stop  3  below) 

8. 

OPTDmi  SOUTH  OF 

B 

n 

(no  stop  until  after  9) 

9. 

raiHAL-DOa  SOUTH OSB 

B 

pi 

(final  print  after  8,  13  or  lk) 

10. 

UHBOUHDKD  SOUTH OF 

S  IS  _ 

A 

pi 

(Terminal  3) 

11. 

RIGHT  HARD  8XDI  QPIR 

S  IS  _ 

0 

pi 

(one  of  2  possible  PIP  termina¬ 
tions,  similar  to  10) 

12. 

MATRIX  8HBULAR 

'  18  _ 

A  (blank)  (Inversion  or  arbitrary  trans 

13. 

BAflr  IFVHRTKL) 

B 

m 

(proper  termination  for  Inver¬ 
sion,  no  stop  until  after  9) 

ftU. 

THECA  AT  MAXDCM 

B 

(the  other  possible  PIP  tormina 

ticn,  similar  to  8;  no  stop 
until  after  9) 


gKggl  3WITCHP 
(If  on) 

1.  Farce  end  -at -stage  and  halt. 

2.  Clslt  the  d^  check  In  Code  5. 

3.  Put  all  print  output  on  tape  6,  Including  any  printed  on-line,  and 
do  a  cycle  print  (to  tape)  every  Iteration. 

k.  On-line  cycle  print. 

3.  Recognita  "  curtain"  narke  In  . 

6.  a.  Vhile  loading  or  after  etops  2  or  3»  load  value  of  f  'from  card. 

b.  Vhlle  running:  perform  a  check  print  at  end  at  Iteration  and 
halt.  (This  le  Independent  of  cycle  print  and  follow*  It.) 

c.  Vhen  using  Code  19  (automatic  pick-up),  load  nev  JO  and  halt. 
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fh ere  are  a any  er^or  rtope  built  into  the  code*.  In  order  to  facilitate 


de-bugglng  of  problems,  they  hare  been  left  separate  and  distinct  and  identi¬ 


fied  vlth  a  auixtr .  Due  to  the  vay  the  codes  are  organized,  it  *4#  not  feas¬ 
ible  to  Identify  these  stops  by  the  B38  location.  In  any  event,  t£is  would 
not  pemlt  easy  reference.  Instead,  the  stops  are  all  "Balt  and  Proceed" 


Orders  vlth  the  stop  mst>er  in  the  address.  Bence  when  the  nachlne  stops, 


the  identifying  ember  appears  in  the  address  part  of  the  Storage  Register. 


All  references  to  these  stops  are  in  octal. 


Brror  Stops  In  Data  Ass 


4444  Tape  6  (input)  will  not  read  properly  after  5  trlee. 
r>252  Zero  everpunch  on  non- zero  digit  in  aj. 

6666  |  a*  j  =£  8192 . 

7700  i>a  for  a*,. 

7770  i>n  for  b1  or  c1  . 

7777  i>*  for  basis  heading. 

10421  Tape  5  (just  written)  will  not  read  properly  after  5  tries. 

This  Is  only  checked  when  the  Matrix  is  punched  in  cards. 


ftwse  are  syetenlsed  as  foliar*. 

Less  than  100  -  Legitimate  stops  in  the  KIR 
100  series  -  Xrror  stop*  in  the  ICR 

200  series  -  "  "  during  pricing  operation*  (Codes  4,  5) 

300  series  -  * '  "in  Code  7  (transformation  of  column) 

400  series  •  "  "in  Cade  9  (making  .change  of  basis,  etc.) 

500  aeries  -  "  "in  Code  11  (recomputing  right  hand  side) 

600  series  -  "  "  during  loading. 

700  series  -  "  "  from  tape  check  failures. 


The  emulate  list  la  In  the  following  table.  The  abbreviation  e-o-s  is  used 
for  end -of-s tags  and  "S«x/S7ARTM  for  "  Put  Switch  x  dosra  before  hitting 
gEAFT  button  I' 
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IS  6 
i  I  S 

MO  O 


S  S  S  S  I  S  S  S  S  S  S  S-  1 


I  1 1  i  e 

o  o  o  o  X 


ii i i i t I : 

OOO  9  l_>  9  R  A 


III 

ill 


I 1 

X  X 


8  8 

o  o 

3  3 

SO  B 
4->  4J 

§  s 

rH  H 

3  3 


.  1 

II 

a : 

•  ►-» 
•*  • 

II 

rH  <3 

€  S 


■  l  * 

iK' 

t  %  1 

!  1 : 
o.  rH  A 

v  5  A 
|58 
Mo 

A  ■ 

X  v. 

V  f 

s  » 

!  1 1 
«  1  « 

III 

5b! 
a  fc  • 

B  -•  rH 

i  a 
^  *  * 


O  rH 

I 

■  10 

**  • 


x  «o 


,  •  8  8 

■g  r 

8  M 

^  1 1 
s  £  - 

1 1 5 

2 1 1 

p  b 

0  •• 

B  rH  w 

Ti  ,0 

p.  3  i 

%  »• ! 

.  8  « 

ff 

♦»  3  -H 

P  rH 

9  0S 


w 

S  I 
k  5  „ 

x  -*  “cT 

tJ  rH  I 

X  H  *a 
«»  w  VO 

I  i  t 

;h 

X  A 
X  • 

til 

j » i 

15  £  c 

rH  jB  +> 


1  i  I 

^  ’8  +* 

2  9  i 

f  8  I 

t :  s 

X  I* 

§  s  ? 

a  *  a 

g  I 1 

f  |2 

:  i § 


*.  8 

ti  «H 


4*  V.«J«©03-H|iO 

a  a « *  £  u  3  s  J I 


1 1 8  i  l 


o  o 


! ;  1 1 
^ « I  i 


8  . 
8  § 

ll 

b  ■* 


•  •  b  t; 

Ml'? 


X  o 

■o  >o 

5  tf 

-fH  <H 

ft  I 


1  i 

3  I 

i* 

H  • 

IS 


e  ►> 

o  ♦> 


5  *8  8 

a  —•  HJ 

h  Vt  o 

•p  *H  « 

S  I  * 


p.k6 


I  5 


|  b  «  fc  * 

I  ■a  I  I  s 

O  fc  H 

*  i  $  a  f 

li  H  X  -rH 


"III 


I 


J  %-t 

£  8 


-885 

sl  U 

o  a  8 

MM 

M  ^  i- 

K  o  to 

jiff 

■P  ®  I 

1  a  1 1 

6  s  -* 


S  S 

S  -n 


O  rH  CM 


•xt  [T\  vo  t-  o  ir»  t— 


u  & 


>o, 

O  O 
rH  rH 


? 

IS  S  S 

ms 


J* 

-* 

IT\ 

«Tv 

t- 

t— 
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0 

0 

0 

0 

# 

0 

0  0  0 

■8 

0 

8 

0 

U 

CJ 

8 

0 

8 

0 

8 

0 

888 

CJ  0  0 
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♦* 

1 

0! 


5 
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V. 
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8 
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a 
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3 
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a 

« 

► 

H 

ta 
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§ 


8 

6 

3 

3 

8 

B 
3 
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3 
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I 
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% 
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Ti 
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3 
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04 


o 

a 

8 

rH 

S 
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O 

ft 


c 

o 

< 

a 

I 


8 

o 

© 
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8 
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c 

0 

► 

0 

2 
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3 

V. 


<5 
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0 

*8 

KJ 

a 

M 


*3 

6 


8 

•H 
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0 

o 


p. 
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8. 

I 


Pi 

o 


r-  o 

»H 

»rv 
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ic 


8 

O 

rH 

W\ 


I 

& 


8 

rl 

♦» 

0 

O 


'3. 


I  3 

I  H 

!  8 
I  0 
I 
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3 
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C 
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8: 
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0 


«-  - 
8*  * 
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v 

3 


8 

V 
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1 


K 

•r-l 
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a> 

1 


<D 


•H  04 
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K\  U-N 

2  i 


N. 

+ 


8 
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p* 
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CD 

I 
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rH 

-0 
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$ 

rH 

i 
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1  A  check  run  error  in  the  loadii*  of  ZJFBASZ  vlll  cause  &  helt  at  7776  r  'tal  with  no  epeclal  indlnatiou 
Hit  SZABT  to  continue  loading.  Ignoring  error.  The  aaount  of  error  if  la  the  Accumulator  on  all  check 
era  rtepe.  Tfclng  Code  19  to  load  a  new  ICR  (vpeclal  ewitch-orer  procedure),  a  check  eta  error  In  the 
carda  caueee  a  atop  with  no  epecial  Indication.  Should  this  occur.  It  veuld  easily  he  overlooked. 


The  remainder  of  the  tape  check  err  or  steps  occur  la  Cods  15*  little  START 
after  these  stops  causes  the  progress  to  continue  as  though  there  tad  been  no  error,  at 
the  operator's  risk. 

762  Rev  record  just  vrltten  on  tape  2  vill  not  read  correctly. 

765  Rsvly  vrltten  tape  5  nay  be  no  good  (tape  check  on  tape  k) 

. . k . .  "  "  "  5) 

765  7fe  and  765  combined. 

766  762  and  l6k  combined. 

There  are  tve  other  error  stop#. 

77777  Automatic  restart  trm  drus  (Code  19)  tried  too  early  during  Inversion 
run.  Reload  deck. 

717777  Cede  19  got  a  check  sun  error  In  restoring  338  fron  dnsss .  This  check  sun 
Is  carried  on  K-,  H-,  and  T-reglons  only.  The  first  7  In  this  error  rnafcar 
is  In  the  tag  field  of  the  Storage  Register. 


J 


RtBOTOO  TOflB 

It  la  Impossible  to  giro  any  realistic  tins  astlnates  on  a  Job.  Tim  per 
Iteration  goea  up  with  a,  with  a  and  vlth  the  mnfcer  of  iterations  alroady  par* 
for— d .  Klaborate  formula  hare  boon  vorked  out  for  thla  but  thoy  art  not  worth 
writing  down  although  they  hart  provided  a  cbm  cooperative  information.  lor  la 

paat  experience  a  reliable  guide.  Aa  lor  era  Ion  run  on  a  195  order  system  (a  •  195) 

t  . 

baa  been  observed  clipping  along  •t  1%  aeca.  per  Iteration.  A  51  order  aywtm 
took  2#  nine,  to  reach  an  opt  Inal  oolutloa.  Clearly  three  numbers  da  not  bear  a 
elaple  relation  to  one  another.  One  of  the  critical  factors  le  the  denalty  at  the 
natrix.  The  following  atateaanta  can  be  node,  hewer  ex . 

1.  The  inversion  NBR  la  the  faeteet  In  operation,  at  leaet  twice  aa  faet  aa 
the  COM  1CF  far  a  given  problem. 

2.  The  COMf  NCR  perform  elaplex  cycle#  faster  than  any  other  known  code 
considering  the  alee  of  problem  allowed,  the  accuracy  maintained  and 
the  flexibility  provided  for. 

3.  PLP  is  invariably  slower  after  a  given  mnfeer  of  iteration#  than  the 
COMM  NCR.  It  is  probably  vlae  to  re -Invert  the  optimal  basis  dbtaimd 
by  the  COMM  NCR  before  beginning  extended  calculation#  vlth  IIP. 


